AD-A264  416 

I  (Mini  llifcn.'i  V 


•ii'lfillilillil  |!| 


Vivacity  Calculation,  Rolled  Ball,  and 
Rosette  Grain  Form  Function  and  Editing 
Capability  Update  to  the  Closed -Chamber 
Data  Analysis  Program  BRLCB 

William  F.  Oberle 


ARL-MR-61  April  1993 


DTJC  1 

i'^l-ECTEWk,' 

iqgrl  Bf 

%/  B  y 


APPROVED  FOR  PUBIJC  RELEASE;  UISTIUBimON  IS  UNUMITED. 


93-10925 


93  5  14 


NOTICES 


Destroy  this  report  when  it  is  no  longer  needed.  DO  NOT  return  it  to  the  originator. 


Additional  copies  of  this  report  may  be  obtained  from  the  National  Technical  Information 
Service,  U.S.  Department  of  Commerce.  5285  Port  Royal  Road.  Springfield,  VA  22161. 


Tt.e  findings  of  this  report  are  not  to  be  construed  as  an  official  Department  of  the  Army 
position,  unless  so  designated  by  other  authorized  documents. 


The  use  of  trade  names  or  manufacturers’  names  in  this  report  does  not  constitute 
indorsement  of  any  commercial  product. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


PuDlic  '^pooftinq  ouropn  tor  tnis  collection  of  intorm,»Cior  is  ‘'Stim<»teo  to  ^-eriqe  i  t'our  oer  roipoese.  including  tne  time  for  fevi»»'«vinq  mstroctioes.  se-^rcrunq  «>..stinq  ^ourc^s 
>)<«thennq  jna  nMintaminq  the  data  needed,  jnj  comoietinq  jna  reviewinq  the  •  oilection  of  information  Send  comments  regarding  th'S  bcifden  estimate  or  mv  jther  .isoect  Of  tnis 
collection  of  information,  including  suggestions  for  reducing  this  Durden  to  Washington  Heaaouarters  Services,  Directorate  for  information  Ooeratiof's  and  fteoorts.  12 1 S  jetferson 
Oavis  Highway.  Suite  1 204.  Arlington.  vA  22202-4102  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0 1 88).  Washington.  OC  20S03 


1.  AGENCY  USE  ONLY  (Leave  blar\k)  I  2.  REPORT  DATE 

I  Anril  1993 


4.  TITLE  AND  SUBTITLE 


3.  REPORT  TYPE  AND  DATES  COVERED 

Final,  Jun  92  -  Jan  93 


5.  FUNDING  NUMBERS 


Vivacity  Calculation,  Rolled  Ball,  and  Rosette  Grain  Poim  Function  and  Editing 
Capability  Update  to  the  Closed-Chamber  Data  Analysis  Program  BRLCB 


6.  AUrHOR(S) 


DA3 11880 

IF2Z9W 

9XDGS3 


William  F.  Oberlc 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  AOORESS(ES) 

U.S.  Army  Research  Laboratory 

ATTN;  AlvlSRL-WT-PA 

Aberdeen  Proving  Ground,  MD  21005-5066 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


t  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AOORESSlES) 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRL-OP-CI-B  (Tech  Lib) 

Aberdeen  Proving  Ground,  MD  21005-5066 


10.  SPONSORING/ MONITORING 
AGENCY  REPORT  NUMBER 


ARL-MR-61 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


12b.  DISiRIBUTION  CODE 


13,  ABSTRACT  (Maximum  ilOO  words) 

Updates  to  the  closed-chamber  data  analysis  program  BRLCB  are  documented.  The  updates  include  vivacity 
calculations,  grain  geometry  form  functions  for  rosette  and  rolled  ball  propellants,  aiid  improved  editing  capabiiities. 


14.  SUBJECT  TERMS 


closed-bomb,  closed-chamber,  *  uming  rate,  propellants,  rolled  ball  propellant 


IS.  NUMBER  OF  PAGES 

43 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  MB.  SECURITY  CLASSIFICATION  119.  SECURITY  CLASSIFICATION  J  20,  LIMITATION  OF  ABSTRACT 

OF  REPOKf  I  OF  THIS  PAGE  I  OF  ABSTRACT  I 


UNCLASSIFIED 


MSN  7'.L0.01.230-S‘,00 


UNCLASSIFIED 


UNCLASSIFIED 


Stand.-ird  Form  29B  ..Rev  2-39) 

Pr^Uftbcd  by  ANSI  Sta  .'I')- '8 


Intentionally  left  blank. 


11 


TABLE  OF  CONTENTS 


1.  INTRODUC'TION  .  1 

2.  VIVACITY  CALCULATION  .  1 

3.  ROLLED  BALL  AND  ROSETl’E  PROPELLANT  GRAIN  FORM  FUNCTION  ....  3 

4.  EDITING  CAPABILITY .  8 

APPENDIX  A;  LISTING  OF  CODE  ADDED  TO  BRLCB  FOR 

VIVACITY  CALCULATION  .  11 

APPENDIX  B:  LISTING  -  PROGRAM  EDITIN.FOR  .  15 

APPENDIX  C:  LISTING  -  PROGRAM  EDOUT.FOR  .  19 

DISTRIBUTION  LIST  .  37 


LIST  OF  FIGURFS 


Figure  Page 

1.  Percent  difference  between  BRLCB  dynamic  vivacity  calculation  and 

direct  calculation  of  the  dynamic  vivacity  from  the  pressure-time  data .  4 

2.  Schematic  of  roiled  ball  grain  geometry .  5 

3.  Percent  difference  between  BRLCB  derived  bum  rates  and  assumed  bum 

rate,  rolled  ball  fonn  function .  6 

4.  Schematic  of  rosette  grain  gcomcti^ .  7 

5.  Percent  difference  between  BRLCB  derived  bum  rates  and  assumed 

bum  rate,  rosette  fonn  function .  9 


LIST  OF  TABLES 


Table  Page 

1.  Format  of  Vivacity  File .  2 

2.  Input  Parameters  Utilized  in  IBHVG2  to  Generate  a  Pressure  History 

Used  to  Validate  the  BRLCB  Vivacity  Calculation .  2 

3.  Input  Parameters  Used  in  IBHVG2  to  Generate  a  Pressure  History 

for  Rolled  Ball  Grain  Geometry  .  5 

4.  Relation  Between  Number  of  Rings  and  Perforations  .  7 

5.  Input  Parameters  Used  in  IBHVG2  to  Generate  a  Pressure  Histi  /  for 

Rosette  Grain  Geometry .  H 

ft.  Example  of  EDITFL  File  .  10 


V 


INTENTIONALLY  LEPT  BLANK. 


VI 


1.  INTRODUCTION 


Although  the  Closed- Chamber  Data  Analysis  Program  BRUCB  (Oberie  ;ind  Kookcr  1993)'  was 
originally  developed  to  support  the  work  of  the  High  Pressure  Combustion  Research  Team  of  the 
Propulsion  and  Flight  Division  of  the  U.S.  Army  Research  Laboratory,  its  use  throughout  the  combustion 
community  has  resulted  in  the  identification  of  additional  features  which  if  incorporated  would  enhance 
the  programs  utility.  The  requested  features  were:  (1)  calculations  to  determine  vivacity,  (2)  inclusion 
of  additional  propellant  grain  form  functions  for  rolled  ball  and  rosette  geometries,  and  (3)  editing 
capabilities  when  dealing  with  multi-layered  propellant  grains  (specifically,  the  ability  to  edit  the  input 
thcrmochcmical  information).  The  objective  of  this  report  is  to  document  the  inclusion  of  these  options 
into  the  BRLCB  program  and  provide  the  user  with  details  on  their  u.se. 

2.  VIVACITY  CALCULATION 

In  addition  to  tlie.  measurement  of  a  propellant’s  bum  rale,  a  compari.son  of  the  closed-chamber 
pressure-time  curve  with  the  pressure-time  curve  of  a  ".standard"  or  ba.sclinc  propellant  is  often  u.sed  to 
characterize  a  propellant’s  combustion  behavior.  Some  of  the  more  common  measuretncnt.s/comparisons 
arc  dynamic  vivacity,  regular  vivacity,  relative  force,  and  relative  quickness.  Generally,  these  comparisons 
or  measurements  involve  the  time  derivative  of  the  prc.s.sure  (dP/dt),  maximum  pressure  (Pmax),  and 
instantaneous  prc.ssurc  (P).  For  example,  dynamic  vivacity  is  defined  by 

Dynamic  Vivacity  = _ -  .  (1) 

(P  *  Pm.'x) 

In  order  to  facilitate  die  calculation  of  the.se  quantities,  an  additional  sub-option  has  been  added  to 
BRLCB.  This  sub-option  is  provided  to  the  u.ser  Just  prior  to  exiting  Option  5,  Smooth  Prc.ssurc/Time 
Data,  of  die  Main  Menu.  Before  exiting  Option  5,  the  u,ser  is  offered  the  option  of  perlonning  the 
vivacity  calculation.  If  the  vivacity  calculation  is  .selected,  the  user  is  a.sked  to  provide  a  llle  name  for 
the  vivacity  infomiation.  The  BRLCB  program  places  no  restrictions  on  the  file  name  fonnai;  however, 
a  file  extension  of  .VIV  is  recommended  to  prevent  duplicating  file  names  automatically  generated  by 
BRLCB.  The  fomiat  of  the  vivacity  tile  is  given  in  Table.  1  and  consists  of  five  columns,  as  shown  in 
the  table.  Although  the  dynamic  vivacity  is  the  only  above  mentioned  quantity  calculated,  :dl  the  other 


Oborle,  W.  I-'.,  aiid  0.  L.  Koukcr.  "BRLCH:  A  Closcd-Clijuiiljcr  Data  Analysis  Program  Pai  t  I  and  It-  Tlicory  and  U.scr's 
Manual."  ARI.-TRA6,  U  S.  Aniiy  Rusuartli  LalKiralory.  Aberduen  IXToving  (iround,  MD,  January  IVd.l. 
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quantities  can  be  computed  using  the  information  supplied  in  the  table,  provided  the  "reference  propellant" 
information  required  for  relative  force  and  quickness  is  known. 


Table  1.  Format  of  Vivacity  File 


Column 

1 

2 

3 

4 

5 

Quantity 

Time 

Pressure 

Prcssurc/Pmax 

dP/dt 

Dynamic  Vivacity 

Units 

s 

MPa 

(••) 

MMPa/s^ 

l/(MPa-s) 

MMPa  -  Million  MPa 


To  validate  the  computer  coding  of  the  vivacity  sub-option,  the  interior  ballistic  code  IBHVG2 
(Anderson  and  Fickle  1987)^  was  used  to  simulate  a  pressure  history  (see  Table  2  for  specific  input  data 
utilized  in  the  simulation). 


Table  2.  Input  Parameters  Utilized  in  IBHVG2  to  Generate  a  Pressure 
History  Used  to  Validate  the  BRLCB  Vivacity  Calculation 


L'hamber  Volume  (cm^) 

300 

Propellant: 

Geometry 

Spherical 

Diameter  (cm) 

0,5 

Mass  (g) 

91 

Impetus  (J/g) 

1,140 

Flame  Temperature  (K) 

3,410 

Molecular  Weight  (-) 

24.87 

Covolume  (cm^/g) 

0.996 

Ratio  of  Specific  Heals  (-) 

1.225 

Density  (g/cm^) 

1,58 

The  prc.'isure  history  generated  by  tBHVG2  was  then  used  in  BRLCB  witlr  the  vivacity  calculation 
option  selected.  Results  from  the  BRLCB  vi^'acity  calculations  were  then  compared  with  results  obtained 
from  directly  computing  the  quantities  listed  in  Table  1  from  Ifie  pressure  history.  The  percentage 


t  Andc/.son,  R.  U.,  unci  K.  13.  Reekie.  "!HHV(i2-  A  User's  (iiiidc."  HRLTR  282V,  U  S.  Army  Rcseureli  LaljorciUny,  Alrcrdccii 
Proving:  (Ironiid.  Ml).  July  19K7. 
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difference  between  tlic  two  calculations  is  shown  in  Figure  1.  As  seen  in  the  figure,  the  maximum  percent 
cnor  is  approximately  0.02%.  Considering  that  dP/dt  was  numerically  calculated  for  both  the  BRLCB 
and  the  direct  calculation  the  results  arc  consider  virtually  identical.  Thus,  the  vivacity  routine  included 
in  BRLCB  is  correctly  performing  the  desired  calculations.  Appendix  A  contains  the  code  listing  for  the 
vivacity  calculation  added  to  BRLCB. 

3.  ROLLED  BALL  AND  ROSETTE  PROPELLANT  GRAIN  FORM  FUNCTION 

Two  additional  propellant  grain  form  functions  have  been  added  to  BRLCB — die  rolled  or  squashed 
ball  and  a  general  ro.settc  form  function. 

The  geometry  of  the  rolled  ball  is  shown  in  Figure  2.  Required  dimensions  are  the  overall  grain 
ditimeter,  D  in  Figure  2,  and  the  thickness  between  the  Hat  surfaces  of  the  grain,  T  in  Figure  2.  The 
equations  utilized  to  compute  grain  volume  and  surface  area  arc  those  used  in  the  interior  ballistic  code 
1BHVG2. 

As  described  in  die  User’s  Manual  (Oberlc  and  Kooker  1993),^  to  verify  the  :  oding  for  a  grain  fomi 
function,  1BHVG2  is  used  to  generate  a  pressure  history  for  a  simulated  closed-chamber  firing.  Input 
values  utilized  in  IBHVG2  for  the  roiled  ball  grain  geometry  arc  given  in  Table  3. 

The  percent  difference  between  the  BRLCB  derived  bum  rate  and  the  assumed  bum  rate  used  in 
IBHVG2  (.see  TabL  3)  is  shown  in  Figure  3.  Due  to  die  decimal  accuracy  of  the  input  data  from 
IBHVG2  and  roundoff/truncation,  errors  as.sociated  with  the  numerical  routines  utilized  in  die  program 
the  maximum  percent  difference  of  less  than  2%  is  consider  insignificant  (see  User’s  Manual  lOberle  and 
Kooker  1993)'  for  details).  Thus,  the  coding  for  the  rolled  ball  propeliant  is  correct. 

The  geometry  of  the  rosette  form  funetion  is  shown  in  Figure  4.  The  required  dimensions  are  outer 
diameter;  lengdi;  perfonition  tliameter;  inner,  middle,  and  outer  web.  ‘  lowever,  ior  the  equations  used  lo 
detennine  volume  and  surlacc  area  (identical  to  tho.se  used  in  1BHVG2),  the  inner  and  middle  webs  must 
be  equal.  In  addition,  the  form  function  is  generalized  for  7,  19,  37,  etc.,  perforations  which  arc 
detemiincd  by  entering  the  number  of  perforation  rings.  The  lelation  between  rings  and  the  most  common 
number  of  perforations  is  given  in  Table  4.  The  generalized  rclation.siiip  (n  is  number  of  rings  and  P 

t  Obcric,  W.  I'.,  :ui(l  D.  E.  Kooker.  "HRlXlIt:  A  t'lo.<:cd-('li:i,nibi-.r  Data  Analy.sis  Pro^^rani  I’arl  1  and  II  -  Tlieory  and  U.scr’.s 
Manila!."  ARL.TR  16.  Ib.S.  Army  Rivsearcli  Lahoialory.  Aberdeen  I’rovinf;  (bound.  Ml),  January  1191. 
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0.02 


Figure  1.  Percent  difference  between  BRLCB  dynamic  vivacity  calculatiun  and  direct  calculation  of  the  dynarnic  vivacity  frgni-tiig 
pressure-time  data. 


D:  Overall  Grain  Diameter 


Tabic  3 


T:  Grain  Thickness 


Figure  2.  Schematic  of  rolled  ball  izrain  geometry. 


Input  Parameters  Used  in  IBHVG2  to  Generate  a  Pressure  History  for  Rolled  Ball 
Grain  Geometry 


Chamber  Volume  (cm^) 

3(X) 

Propellant: 

Geometry 

Rolled  Ball 

Overall  Diameter  (cm) 

0.0 1 

Thickness  (cm) 

0.(K)5 

Mass  (g) 

91 

Impetus  (J/g) 

1,140 

Flame  Temperature  (K) 

3,410 

Molecular  Weight  (g/g-mol) 

24.87 

Covolumc  (cmVg) 

0.99h 

Ratio  of  Specific  Heals  (-) 

1.225 

Density  (g/cm^) 

1.58 

Bum  Rale  Law 

Coefficient  (cm/s-MPa") 

O.lh 

Exponent  (-) 

0.908 
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Figure  3.  Percent  difference  beiween  BRLCB  derived  bum  rates  and  assur;  :(i  bum  rate,  rolled  ball  form  function 


number  of  perforations)  is 

P  =  3 »n^  +  3 *n  +  !  . 

A  similar  procedure  to  that  iscd  to  validate  the  coding  for  the  roiled  ball  form  function  was  used  for  the 
rosette  form  function.  Input  values  for  the  IBHVG2  simulation  arc  given  in  Table  5.  The  percent 
differences  between  derived  (BRLCB)  and  assumed  bum  rates  arc  shown  in  Figure  5.  Although  the 
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maximum  percent  difference  is  slightly  higher  than  for  the  rolled  ball  the  results  arc  felt  to  be  acceptable. 
Thus,  the  rosette  fonn  function  appears  to  be  correctly  coded. 

Table  5.  Input  Parameters  Used  in  IBHVG2  to  Generate  a  Pressure  History  for  Rosette 
Grain  Geometry 


Chamber  Volume  (cm^) 

300 

Propellant: 

Geometry 

Rosette 

Diameter  (cm) 

2.9 

Length  (cm) 

3.5 

Perf  Diameter  (cm) 

0.3 

Inner  &  Middle  Webs  (cm) 

0.2 

Outer  Web  (cm) 

0.3 

Mass  (g) 

91 

Impetus  (J/g) 

1,140 

Flame  Temperature  (K) 

3,410 

Molecular  Weight  (-) 

24,87 

Covolumc  (cm^/g) 

0.996 

Ratio  of  Specific  Heats  (-) 

1.225 

Density  (g/cm^) 

1.58 

Bum  Rate  Law 

Coefficient  (cn/s-MPa”) 

0.16 

Exponent  (-) 

0.908 

4.  EDITING  CAPABILITY 

For  multi-layered  propellant  grains,  12  different  values  (propellant  density  and  ihcrmochemical 
infonnation)  arc  required  by  the  code  for  each  layer.  Such  a  process  is  extremely  siiscentihle  to  error 
(typing,  transposing  digits,  etc.).  Unfortunately,  to  correct  an  error,  the  original  version  of  BRLCB 
required  the  user  to  re-enter  all  12  values  for  each  layer,  again  creating  a  situation  which  is  extremely 
error  prone.  To  alleviate  this  difficulty,  an  additional  option  (Option  0)  has  been  incoiporatcd  into  the 
BRLCB  program.  This  option  will  alUw  the  user  to  edit  the  physical  properties  (density  and 
thermochemistry)  of  the  propellant. 
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igure  5.  Percent  difference  between  BRLCB  derived  bum  rates  and  assumed  bum  rate,  rosette  form  hinction. 


The  option  is  accessed  from  the  Main  Menu  by  typing  a  zero.  Tiic  user  is  then  requested  to  supply 
the  name  of  the  Master  Information  File  (.MAS  extension)  containing  the  propellant  data  to  be  conected. 
After  entering  the  file  name,  the  user  will  be  allowed  to  edit  the  physical  properties  of  each  layer  of  the 
propellant  material.  In  operation,  the  program  reads  die  indicated  Master  Information  File,  creates  a  file 
named  EDITFL  with  the  format  shown  in  Table  6  (example  is  a  three-layer  grain)  and  then  invokes  the 
EMACS  editor  JOVE  (program  EDITIN.FOR,  .see  Appendix  B).  After  the  user  completes  editing  (.lOVE 
is  exited  by  typing  ESC-Z,  escapie  key  followed  by  the  letter  /.),  die  Master  Information  File  is  rewritten 
(program  EDOUT.FOR,  .see  Appendix  C).  THE  NEW  INFORMATION  WILL  OVERWRITE  THE 
ORIGINAL  MASTER  INFORMATION  INDICATED  BY  THE  USER.  THIS  OPFION  SIMPLY  EDITS 
AN  EXISTING  MASTER  INFORMATION  FILE  AND  DOES  NOT  CREATE  A  NEW  FILE. 


Table  6.  Example  of  EDITFL  File 


Number  of  Layers:  3 

LAYER  #;  1 

Property 

Beginning  Layer 

End  of  Layer 

Impetus  (J/g) 

1,100.00 

1,200.00 

Flame  Temp.  (K) 

3,100.00 

3,2(X).(X) 

Density  (g/cm^) 

1.500000 

l.bOOIXX) 

Mol,  Wt. 

24.(XX)000 

25.000000 

Covolumc  (cm^/g) 

l.OOOO(K) 

14.100000 

Gamma 

1.2000(X) 

1.2100a) 

LAYER  #:  2 

Property 

Beginning  Layer 

End  ol  Layer 

Impetus  (J/g) 

L2(X).0() 

1 ,250.00 

Flame  Temp.  (K) 

3.2{X).(X) 

3,300.(X) 

Density  (g/cm^) 

L4(.XXXX) 

1.450000 

Mol.  Wt. 

25.(XXXX)0 

27.(XXXXX) 

Covolumc  (cin-^/g) 

1.2(XXXX) 

1.210000 

Gamma 

1.3(XXXX) 

1.310000 

LAYER  #:  3 

Property 

Beginning  Layer 

End  of  Layer 

Impetus  (J/g) 

1,300.(X) 

1,400.00 

Flame  Temp.  (K) 

3,400.00 

3,450.00 

Density  (g/cm^) 

l.4(XXXX) 

1 .500(X)0 

Mol.  Wt. 

22.00(XXX) 

21.0(XXXX) 

Covolumc  (cm^/g) 

1.3rXMXX) 

1.4(XX)00 

Gamma 

L22(XX)0 

1.250{XX) 
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APPENDIX  A; 


LISTING  OF  CODE  ADDED  TO  BRLCB  FOR  VIVACITY  CALCULATION 


II 


lNT;.NTIONALLY  LEFT  BLANK. 
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(^★★★★★★★:^A***-Ar*******  **•*+•**•**  ■*,■•*:■*••*:*********■****■** -A- -A-****** -A- *3lr**Tlr'Ar** 

^•**********AA ** ***  VIVACITY  CALCULATION  CAN  BE  ADDED  HERE  ******* 

(J  A.***  ****A*Ar**A:*:**:***:  A- Ar-AAA'*******'****'*-***-*^-*'**********  ********* 

353  CALL  CLEAR 

WRITE ( *, 2781) 

2781  FORMAT {///,'  Do  you  wish  to  compute  vivacity  information?', 
1/,'  (l=Ye3,  2=No) ' ) 

WRITE (*,*)' Enter  your  choice.' 

READ (*, *) IKLB 
IF  (IKLB  .EQ.  1)  THEN 

WRITE (*,*)' Enter  a  file  name  for  the  vivacity  information.' 
READ (*, 2782) VIVFIL 
27  82  FORM.AT(A2  0) 

OPEN (UNIT=8, FILE=VIVFIL) 

WRITE (8, 2783) 

2783  FORMAT ('  Time  (sec)  Pressure  (MPa)  P/PMAX  DP/DT 

1  Dynamic  Vivacity' , 

2/,'  DP/DT  has  units  of  million  MDa/sec', 

3/,'  Dynamic.  Vivacity  has  units  1/ (MPa*sec )  '  , 

4/^  /  **A**A,V*AAAA*AA***AA************AAAAAAAA*AAAAAAAA  +  AaAAAA^  ^ 

5//) 

END  IF 

OPEN (UNIT-3,  FILE  =  Al(14),  STATUS  =  ' NEW' , ERR=8 6 ) 

REWIND (3) 

IST=1 

IF  (NP  .GT.  999)  THEN 

IST=HP-990 

ENDIF 

DO  60  I  -1ST,  NP 
TX=(I-IST) *A3 (42) /lOOO . 

C  TX"T(I) 

IF  (ITYPE  .EQ.  5)  THEN 

WRITE  (3,  *)TX,PR(I)  ,PDOT(I)  ,EF.(I) 

ELSE 

WRITE (3,*)  TX,PR(I) ,PDOT (I) 

ENDIF 

IF  (IKLB  .EQ.  1)  THEM 
POPMAX“PR(I) /A3 (43) 

VIVY=1 .E6*PDOT (I) / (PR(I) *A3(43) ) 

WRITE(8,2784)TX,PR(I) , POPMAX, PDOT ( I ) ,VTVY 

2784  FORPJAT  (5Fi6  .  6) 

ENDIF 

60  CONTINUE 

CLOSE(UNJT-3) 

CLOSE (UNIT-8) 
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Intentionally  left  blank. 
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APPENDIX  B: 

LISTING  -  PROGRAM  EDmN.FOR 


15 


Intentionally  left  blank. 


!6 


PROGRAM  EDITIN 


CHARACTER* 

20 

FNAMEl 

,  A2  (6) *80,  A1  (20),  EDITFL 

DIMENSION 

A3 

(100) , 

P  (11,  15,  5) 

WRITE 

{*, 

*) 

f  Till s 

option  only  appld.es  to  editing 

the' 

WRITE 

(*, 

*) 

'thermochemistry  of  a  multi-layered  grain.' 

WRITE 

(*, 

*) 

' Other 

editing  can  be  achieved  by  selecting' 

WRITE 

(*, 

*) 

'Option  1,  Create  Master  Information 

File .  ' 

WRITE 

(*, 

*  ) 

1000  CONTINUE 

WRITE 

(*, 

*) 

' Enter 

the  file  name  of  the  maste.r  fi 

le  ' 

WRITE 

(*, 

*) 

' which 

is  to  be  edited.  ' 

WRITE 

{*, 

*) 

' include  drive  and  file  extension  if 

necessary 

IFGZ  --= 

1 

READ  {*',  50  00)  FNAMEl 
5000  FORMAT {A20) 

OPEN  (UNIT  =  9,  FILE  =  FNAMEl,  STATUS  =  'OLD',  ERR  =  1010) 

REWIND  (UNIT  =  9) 

GO  TO  1020 
1010  CONTINUE 

WRITE  (*,  6000) 

6000  FORMAT (2 OX, ' The  indicated  master  file  does  not  exist,  try  again.' 

1) 

PAUSE 

GO  TO  1000 


(^•ie-k-kir-k'ieie’kif'k-irir-kic-kii:ir-k-ki(irk-ki(ir'ki('kie’k-k-kir-k-k-k-k-k-kk-kic‘kA-k'A^^if-k-kic-k-k-k-ki^-k'k’k'A-k-k'^'k'k-k'kk’k-A 

C**********  reading  OLD  MASTER  FILE  TO  BUILD  FROM  ********************** 

1020  CONTINUE 

DO  1030  I  ••=  1,  6 

READ  (9,  5010)  A2  (I) 

1030  CONTINUE 
5010  FORMAT (A80) 

DO  1040  I  =  1,  20 

READ  (9,  5000)  A1  (I) 

1040  CONTINUE 

DO  1050  I  - 
READ  ( 9 , 

1050  CONTINUE 

DO  1060  I  =  1, 

DO  1070  J  = 

DO  1080 
READ 

1080  CONTINUE 

1070  CONTINUE 
1060  CONTINUE 

CLOSE  (UNIT  9) 

Q*+****A*A*  CREATING  FILE  WHICH  WILL  UE  EDITED  ************************* 

QA'AA’A'AA'A  AAA  A  AAA  AAAAaA  AAAAAAAAAAAAAAAAAAA'AAxAAAAAAAAAAAA  A  AAA  AAAAAAA^'AAAAA 


1, 

*) 


100 

A3  (I) 


11 

I, 

K  -• 
(9, 


15 
1, 
*  1 


(I,  J,  K) 


OPEN  (UNIT  =  9,  FILE  =  'EDITFL') 

REWIND  (UNIT  =  9) 

QAAAAAAAAAAAAAAAAArAAAAAAAAAAAA*AAAAAAAAA*AAA*AAAAAAAAAAAAAAAAA'AAAAAAAAAA' 

Q  *  ★  ★  A  A  A  A  A- A  A  A  A  A  A  NOW  INFORMTAION  FOR  EACH  LAYER  *  *  ^- •*■  *  ^  •*  *  *^*  *  *•  * ’^ ’^  •’^  ■^  a- a  a  a  a  a  a  a 

(^AAA‘AA;AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA'AAAA 
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NUMLAY  =  INT  (A3  (4)  +  .5) 

WRITE  (9,  6010)  FNAMEl,  NUMLAY 
6010  FORMAT  {A20 ,  T30 Number  of  Layers  ,T47, 12) 

DO  1090  1=1,  NUMLAY 
WRITE  (9,  *) 

WRITE  (9,  6020)  I 

6020  FORMAT (' LAYER  #:',I2) 

WRITE  (9,  6030) 

6030  FORMAT (' Properity' T20,  '  Beginning  Layer' , T4  0, ' End  of  Layer') 
WRITE  (9,  6040)  P  (2,  I,  1),  P  (2,  I,  2) 

6040  FORMAT {' Impetus  ( J/a ) ' , T20 , FIO . 2 , T40 , FlO  .  2 ) 

WRITE  (9,  6050)  P^{3,  I,  1),  P  (3,  I,  2) 

6050  FORMAT (' Flame  Temp.  (K) ' , T20 , FIO . 2 , T40 , FIO . 2 ) 

WRITE  (9,  6060)  P  (4,  I.  1),  P  (4,  1,  2) 

6060  FORMAT ('Density  ( g/cc ) ' , T20 , FIO . 6, T4 0 , FIO . 6 ) 

WRITE  "(9,  6070)  P  (5,  I,  1),  P  (5,  I,  2) 

6070  FORMAT('Mol.  Wt . ' , T20, FIO . 6, T40, FIO . 6) 

WRITE  (9,  6080)  P  (6,  I,  1),  P  (6,  I,  2) 

6080  FORMAT (' Co volume  ( cc/g ) ' , T20 , FIO . 6,  T40 , FIO . 6 ) 

WRITE  (9,  6090)  P  (7,  I,  1),  P-(7,  I,  2) 

6090  FORMAT (' Gamma' , T20, FIO . 6, T40, FIO . 6) 

1090  CONTINUE 

CLOSE  (UNIT  =  9) 

END 
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APPENDIX  C; 

LISTING  -  PROGRAM  EDOUT.FOR 


19 


INTENTIONALLY  LEFT  BLANu. 


20 


PROGRAM  EDOUT 

COMMON  A3  (100),  P  (11,  15,  5) 

CHARACTER*20  FNAMEl,  hi  (6) *80,  A1  (20),  LINE*80 

ic  "k  if  ic  if  -k  "k  -k  -k  if  -k  -k  -k  it  it  -k  ic  iK  -k  -k  k  -k  -k  -k  -k  "k  k  if  it  -k  "k  -k  -k  -k  -tc  k  "k  "k  -k  -k  i(  -k  -k  -k  -k  -k  -k  -k  -k  if  -k  -k  -k  -k  -k  -k  "k  -k  -k  -k  -k  -k  i(  'k  -k  ^  -k 

C*******************  reading  edited  FILE:  EDITFL  *********************** 

OPEN  (UNIT  =  9,  FILE  =  'EDITFL') 

REWIND  (UNIT  =  9) 

READ  (9,  5000)  FNAMEl,  NUMLAY 
5000  FORMAT (A20,T47, 12) 

OPEN  (UNIT  =  10,  FILE  =  FNAMEl) 

REWIND  (UNIT  =  10) 


(^k  k  k  k  -k  k  k  k  k  k  k  k  if  *.  k  It  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  -k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

Qkkkkkkkkkkkk*kkkk  reading  master  file  to  rebuild  ickkkkkkkkkkirkkkic*****  tr 

Qkkk-kkkkkkkkkkirkMkif-kk'*rkkkkkkJfkkkkkkkkkkkkkkkkkkkkkkkTfi(kkkiekk***kkk*kk*kk 


1000  CONTINUE 

DO  1010  I  =  1,  6 

READ  (10,  5010)  A2  (I) 

1010  CONTINUE 
5010  FORMAT (A80) 

DO  1020  1=1,  20 

READ  (10,  5020)  A1  (I) 

1020  CONTINUE 
5020  FORMAT (A20) 

DO  1030  I  =  1,  100 

READ  (10,  *)  A3  (I) 

1030  CONTINUE 

DO  1040  1=1,  11 
DO  1050  J  =  1,  15 
DO  1060  K  =  1,  5 

READ  (10,  *)  P  (I,  J,  K) 

1060  CONTINUE 

1050  CONTINUE 
1040  CONTINUE 

CLOSE  (UNIT  =  10) 

(]*********************************************************************** 
C******  NOW  INFORMTAION  FOR  EACH  LAYER  IS  READ  FROM  FILE**************** 


DO  1070 

I  ••= 

,  NUMLAY 

READ 

(9, 

5030)  LINE 

5030 

FORMAT 

(A80) 

READ 

(9, 

5030) 

READ 

(9, 

5030) 

READ 

(9, 

5040)  P  (2,  I,  1), 

P 

(2, 

I, 

2) 

5040 

FORMAT 

(T20, 

F10.2,T40,FI0.2) 

READ 

1  Q 
\  r 

5050)  P  (2,,  I,  1), 

P 

(3, 

1, 

2) 

5050 

FORMAT 

(T20, 

FIO .2,T40,F10.2) 

READ 

(9, 

5060)  P  (4,  I,  1), 

P 

(4, 

I, 

2) 

5060 

FORMAT 

(T20, 

F10.6,T40,FL0.6) 

READ 

(9, 

5070)  P  (5,  I,  I), 

P 

(5, 

I, 

2) 

50  70 

FORMAT 

(T20, 

FI0.6,T40,FIO,6) 

READ 

(9, 

5080)  P  (6,  I,  1), 

P 

(6, 

I, 

2) 

5080 

FORMAT 

(T20, 

F10.6,T40,F10.6) 

READ 

(9, 

5090)  P  (7,  I,  1), 

P 

(2, 

I, 

2) 
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5090  FORMAT (T20,F10 .6,T40,F10 .6) 


p 

(9, 

I,  1)  =  1.98717/P 

(5,  I, 

1) 

p 

(9, 

I,  2)  =  1.98717/P 

(5,  I, 

2) 

p 

(10, 

I,  1)  =  P  (9,  I, 

1)  /  (P 

(7, 

I, 

1) 

p 

(10, 

I,  2)  =  P  (9,  I, 

2)  /  (P 

(7, 

I, 

2) 

p 

(11, 

I,  1)  =  P  (10,  I, 

1)*P 

(3, 

I, 

1) 

p 

(11, 

I,  2)  =  P  (10,  I, 

2)*P 

(3, 

I, 

2) 

1070  CONTINUE 

CLOSE  (UNIT  =  9) 


1.) 

1.) 


Q-k  -k  -k  -k  -k  -k  -kie  i 


:-k'k'k-kk’k'k’k'ki 


INTEGRALS  ARE  UPDATED 


C** 


kkkkkk-kk-kkki^-kk-k-kk-k'k'kir’kk'kk'k-k-k-k'k-k-k-k-kkk-k'kk-k-k'k-kkkk'k'k'kkk-kk-kk'k'k'kk'k-kkkk'kkk'kk 


CALL  COMMASS 


(^•kkkkkkkk-kkkkkk  kkk-kkk'kkk'kir'k-ir-k'k'k'A'kkick'k'k'k'kirk'kk-k'k'k'k’k’k-k’kickkkkk'kk’k'kkkk'k'k'k'k'kkk 

(2*  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  [^STER  FILE  IS  REWRITTEN  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 


OPEN  (UNIT  =  I,  FILE  =  FNAMEl) 

1030  CONTINUE 

DO  1090  I  =  1,  6 

WRITE  (7,  6000)  A2  (1) 

6000  FORMAT (A80) 

1090  CONTINUE 

DO  1100  :  =  1,  20 

WRITE  (7,  6010)  A1  (I) 

6010  FORMAT (A20) 

1100  CONTINUE 

DO  1110  I  -  1,  100 

WRITE  (7,  *)  A3  (I) 

1110  CONTINUE 

DO  1120  1=1,  11 
DO  1130  J  =  1,  15 
DO  1140  K  -  1,  5 

WRITE  (7,  *)  P  (1,  J,  K) 

1140  CONTINUE 

1130  CONTINUE 
1120  CONTINTIE 

CLOSE  (UNIT  =  7) 

END 

Qk  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

c****  **  **************  *su3R0UTINE  COMMASS*********  *  ********************** 

Q*********************^ ************************************  ************* 


SUBROUTINE  COMMASS 

Q*********** **********************************************  ************** 

C  Version  3.0,  January  1992 
C 

C  THIS  SUBROUTINE  WILL  COMPUTE  THE  INITIAL  MASS  PER  LAYER  AS  WELL  AS 
C  THE  INTEGRALS  OF  EACH  PROPERTY 
C 

Qk  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 


COMMON  A3  (100),  (11,  15,  5) 

IGTYPE  INT  (A3  (45)  +  .5) 

LF  (IGTYPE  .GT.  14)  THEN 
I(;TYPE  =  15 
END  IF 


27. 


:  I 


NL  =  INT  (A3  (4)  +  .5) 

IGQR  =  INT  (A3  (31)  +  .5) 

(J*****A-***'^-iIc**t^**********  **•****************■***★*>•**★* 

C*****************  STORAGE  LOCATION  FOR  INTEGRAL  ARE 

(^i('kii:-k'fc-k-k'f('A'k-k-k-k-k-k-k-ft-kir'k-k'k-k'kiri<^ic'kicitA-fic-k-k-K-k-k1r-k-k-kiriT-k-k-kick'k-*--A 

DO  1000  1=2,  11 
DO  1000  J  =  1,  15 
P  (1,  J,  3)  =  0.0 
1000  CONTINUE 

IF  (IGQR  .EQ.  1)  THEN 
XI  =  0.0 

CALL  FORM!'  (IGTTPE,  ASURl',  VOLUNb,  XI) 

P  (8,  1,  3)  =  VOLUNB 
DO  1010  1=2,  11 

P  (I,  1,  3)  =  P  (I,  1,  1)  *P  (8,  1,  3) 

1010  CONTINUE 

A3  (3)  =  P  (4,  1,  3) 

RETURN 
END  IF 

DO  1020  J  =  1,  NL 

(2 -k  ic  -k  'k  it  "k  ir  it  ic  ir  "k  ic  -A  -k  ic  -k  "k  ir  it  ir  -k  if  -k  ir  ir  -k  -k  if  "k  ir  "k  -k'  "k  "k  ic  "k  "k  ic  ic  -k  ic  ic  it  -k  -k  -k  ic  ic  "k  -k 

C*****************  determine  X  VALUES  FOR  INTEGRATION 

•kkkkkkkkkkkk 


i:itkitkk'kkt\i(kkirkkf'k-k’k 

ZEROED  **  +  *****-v*** 

*■*•★*-***★  VP  V  Tk  tV  *  ★  *  ★ 


kkkkkkkkkkkkkkkkkkk 

kkkkkkkkk 

i  IF  ( J 

.EQ. 

NL) 

THEN 

XS 

=  P 

(1, 

J,  1) 

1  XL 

=  A3 

(1) 

1  ELSE 

i  XS 

=  P 

(1, 

J,  1) 

1  XL 

=  P 

(1, 

J  +  1 

kkkkkkkkkkkkkkk 
•kkkkkkkkkkkiiickkkkk 
kkkkkkkkkkkkkkkkkkk 


1) 

END  IF 

Q^k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

:************  ALWAYS  300  SUBDIVISIONS  ***** 

kkkkkkkkkkkk kkkkkkkkkkkkkkkkkkkkkkkkk 

XDEL  =  XL  -  XS 
XSTEP  -  XDEL/300. 

DO  1030  1=1,  301 

XI  =••  XS  +  (I  -  1)*XSTEP 

(2***********  ****************************  ************ 

CALL  FORMT  (IGTYPE,  ASURF,  VOLUNB,  XI) 

DO  1040  K  “  2,  11 

PDEL  =  P  (K,  J,  2)  -  P  (K,  J,  1) 

^kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk  k  k  k  k 

C**********************  integration  is  PERFORMED  **** 

rkkkkkkkkkkkkkkkkkkkkkkkkkkkk A kkkkkkkkkkkkkkk 

FCN  =  ( fPDEL/XDEL) * (XI  -  XS)  +  P  (K,  J 


kkkkkkkkkkkkkkkkkkk 
kkkkkkkkkkkkkkkkkkk 
k  k  k  k  kkk  ' 


:  k  k k  k k  k k  kk k k  k 


kkkkkkkkkkkkkkkkkkkkkkkkkk 


k  k  kkk  kkkkkkk  kkk  kkk  k 
k  k  k  k  k  k kkkkkk kkkkkkk 
kkkkkkkkkkkkkkkkkkk 


1) ) *ASURF 


ii'  (  (i  .  EQ. 

i) 

.OR. 

(I 

EQ. 

3U1))  THEN 

P  (K,  J, 
ELSE 

3) 

=  P 

(K, 

J,  3) 

FCN 

P  (K,  J, 
END  IF 

3) 

-  P 

(K, 

•J,  3) 

+  2.* FCN 

1040 

1030 


CONTINUE 
CONTINUE 
DO  1050  K  = 


2,  11 
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iiimwBiiiiiiiitiiiiiii 


TmaamwitniiiaMW^^ 


P  (K,  J,  3)  =  P  (K,  J,  3) * (XDEL/600  .  ) 

1050  CONTINUE 
1020  CONTINUE 

Qiti^-kififi^-iciK^'^^ici^iii^icir-k-k'kie-k-k-kkkk-kit-kic-k-k'k-kitic-k-k'k-kic-k-k-k'k'k'k-kic'k'k-k-kkic'k-kk-Aic-kic-k-k-A'k'k-k-k'k 

C:*:*-A:>:**:*****^*:^****:^-***  |yj^SS  QF  SINGLE  GRAIN  *******************-*****  ^' *  * 

QtV****:*:**?^  A-**  A- *A-A***:*r*:A:******  *********************************  A********* 

DO  1060  I  =  1,  NL 

A3  (3)  =A3  (3)  +P  (4,  I,  3) 

1060  CONTINUE 
RETURN 
END 

Q******  *******  A  A  *************************************  ******************* 

C* 

C******Ar***A**A******  SUBROUTINE  FORMT  ********************************* 
C* 

S*********** ***************************************************  ********* 

c* 

C  ICODE:  code  for  type  of  grain 

C  R  :  burn  depth 

C  GL:  unburned  grain  length 

C  D  :  unburned  outer  diameter 

C  PD:  unburned  perforation  diameter 

C  WI,  WM,  WO:  inner,  middle  and  outer  webs  respectively 

C 

C  Output : 

C  SFAREA:  surface  area 

C  FRCSFA:  surface  area/initial  surface  area 

C  VOLUND :  unburned  volume 

C  VOLBRN :  burned  volume 

C  FRCBRN;  burned  volume/ initial  volume 

C  VOLMAO:  unburned  volume  of  outer  layer 

C  VOLMBO:  unburned  volume  of  inner  layer 

C  VOLABR:  burned  volume  of  outer  layer 

C  VOLBBR:  burned  volume  of  inner  layer 

C 

SUBROUTINE  FORMT  (ICODE,  SFAREA,  VOLUNB,  R) 

COMMON  A3  (100),  P  (11,  15,  5) 

DIMENSION  S7  (4),  S19  (3,  4) 

DATA  RT/1 . 732050808/,  PI3/1 . 047197551/,  PI/3.141592654/ 

Q-k  ir  -k  ir  if  -k  -k  -k  *  *  -k  -k  k  *  it  -k  *  -k  -k  ie  A  *  -k  *  -k  *  *  -k  k  ic  *  -k  -k  k  ir  k  ie  ic  *:  ic  k  A’  ic  it  ir  it  ie  it  i>.-  k  *  -k  ic  ic  -k  -k  -k  k  -k  -ic  i(  -k  -k  -k  fc  -k  -k  -k 

C*AArAr*A:AA:**A:A:A:ArA*A-**A:  SET  GRAIN  GEOMETRY  *  *  *  *  *  *  7^  *  A:  *  *  *  *  :A  *  *  *  A  Ac  *  *  *  *  *  *  A:  *  *  *  *  * 

(jAAAAAAAAAAAAAAAAArAAA.AA-AAAAAAArAAAfA-AAA-AArAAAcArAA-AAAAAAA-ArA-AAAAAAAAAfAAAAAArAA:* 


GL 

=  A3 

(7) 

D  = 

A3 

(8) 

PD 

=  A3 

(9) 

WI 

=  A3 

(10) 

WM 

=  A3 

(11) 

WO 

=  A3 

(12) 

(^A’AA’AAArA‘A:A'ArA’A'ArA:A’A’A’A'AA'A:AAArA‘A‘AAtArArArAVtAAfAfArA(ArAfAr.ArArArAr'AA'AArAfAAfA‘;A^Af'frArA’‘AAfA'ArArAr'A*A:A’‘AA 

C********  Set  U  2*  (depth  burned)  and  branch  to  graiti  type  ************ 

(SAfAA'AAArAA-*AAAAAArAAA-A*A:'A’*A’AAAA:AAr*ArA-*A-'AA*ArA*ArAArArA:Ar-AArAr-A*ArA*AA:A:Af**ArArAA:A:A:AfAA-* 

U  =  2 . 0*R 

QAAAAAA’AAArAAAAAA'AAAAAA*A'AAAAAA'AAAfAfAAArA'»VAr*A’*AfArArArAAAAAr*AAAAArAAAcA-A-AAAr'AArArAAA 
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\mmuv6jK.  *-aa4»^ 


£ 


C******  FINDING  NUMBER  OF  RIGNS  FOR  ROSETTE  GRAIN  AND  VALUES  *********** 

{^**-A***Ar*A*****-*:**--Ar***A*************'**A'****:****iki!f***'******-****-**Tt******* 

(^^f'kii(ii(ificiri('kiKi(iiiiir‘ki(icici(ic'kirif-k^it'^'k'k'k-k-iririciric'k1cit‘kiici\‘k'k'k'k'k-!cirif-k'ki(-k'k'kk-k-k-ic-k-kirir-kic-k-kic'kir 

C***********  t’INNDING  NUMBER  OF  RIGNS  FOR  ROSETTE  GRAIN  **  + ^- *********** * 

Q****** ************************ *********************************A ******* 

IF  (ICODE  .EQ.  15)  THEN 

NRGS  -■  INT  (A3  (45)  -  150.  +  .5) 

NIZ  =  6*NRGS*NRGS 
NOZ  =  6*NRGS 
NCZ  -  6 

(^^^★^^•★★★★★Tlr*****^****’'*'****:***^***^  A- *:Arir*********'A:****************'^i:**  ****** 

C*********** A ****  test  to  DETERMINE  IF  METHOD  APPLICABLE  *************** 

Q********  ****************************************************  ************ 

TOPP  =  WO*-  (PD  +  WO) 

BOTT  =  .5*WI*(PD  +  WI) 

IF  (BOTT  .GT.  TOPP)  THEN 
WRITE  (*,  *) 

1  'The  grain  geometry  does  not  satisfy  the  assumption' 

WRITE  (*,  *) 

1  'made  in  determining  the  surface  area  and  unburned' 

WRITE  (*,  *) 

1  'volume  of  a  grain.  Use  any  results  with  caution.' 

END  IF 

DPRIME  =  PD  +  U 

(J  *******************************-****  ************************  ************ 

C*A**************  OUTER  SLIVER  IS  COMPUTED  **********************  >^  ***** * 

AAkiiifi^k-k-k-^ic-i(k-k'k'k'Ak-kfc‘k-kirir'^kir-A'‘kic-k-kicic'kiciric*:ki<i(icirkicicic-Mi(ickit'kic-k‘k'A-k'k‘k-Ai('ki('k-kiri( 

THETA  ••=  2.*ACOS  (AMINl  (1.,  (WI  +  PD) /DPRIME)) 

TOPP  =  WI/2.  +  PD/2. 

BOTT  =  WO  +  PD/2 , 

ALPHA  -  ACOS  (TOPP /BOTT) 

IF  (U  .LE.  WO)  THEN 

AREAl  -  .5*  (PI/2.  -  ALPHA)*- ((PD  -  U  +  2  .  *  WO )  /  2  .  )  *  *  2 
AREA2  =  .5*((WI  +  PD) /2 . ) * (PD/2 .  +  WO) *SIN  (ALPHA) 

AREA3  =  .5*  (PI/2.  -  ALPHA) * (U/2 .) **2 

AREA4  -  .5*  (PI/2.  -  (THETA/2 .))* (DPRIME/2 .) **2 

AREAS  =  .5*((WI  H-  PD) /2  .)*  (DPRIME/2  .)  *SIN  (THETA/2.) 

EOZ  =  2.* (AREAl  +  AREA2  -  AREA3  -  AREA4  -  AREAS) 

XLENl  -  (PI/2.  -  ALPHA)-*  ((PD  -  U  +  2.*WO)/2.) 

XLEN2  (PI/2.  -  ALPHA)  *  (U/2  .  ) 

XLEN3  -=  (PI/2.  -  (THETA/2  .))*  (DPRIME/2  .  ) 

SOZ  -  2.*(GL  -  U)* (XLENl  +  XLEN2  +  XLEN3) 

ELSE 

TOPP  ^  2.*WO*(PD  H-  WO) 

BOTT  =  PD  +  SORT  ((2.*WO  -  WI)*(2.*PD  +■  2  .  *WO  +  WI)) 

IF  (U  .LE.  (TOPP/BOTT) )  THEN 

TOPP  -  PD*DPRIME  *  2.*WO*(Pl)  +  WO) 

BOTT  =  DPRIME* (PD  )  2 . *WO) 

BETA  =  ACOS  (TOPP/BOTT) 

TOPP  =  2.*WO*(PD  H  WO)  -  (PD*U) 

BOTT  =  U*(PD  +  2.*WO) 

GAMMA  =  ACOS  (TOPP/BOTT) 

AREAl  -  .5*  (PD/2.  H-  WO)  *  (DPRIME/2  .)  * 

1  SIN  (ALPHA  -  THET.V2.) 
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AREA2 

=  .5*(PD/2. 

AREA3 

=  .5* (ALPHA 

AREA4 

=  .5*(PI/2. 

EOZ  = 

2 . * (AREAl  - 

XLENl 

=  (PI/2.  - 

XLEN2 

( ALP  HA  - 

SOZ  = 

2 . * (GL  -  U) 

ELSE 

EOZ  = 

0.0 

SOZ  = 

0.0 

END  IF 

END  IF 


+  WO) * (DPRIME/2 . ) *SIN  (BETA) 

-  BETA  -  THETA/2 .)* (DPRIME/2 ,) **2 

-  ALPHA  -  GAMMA) * (U/2. ) **2 
AREA2  -  AREA3  -  AREA4) 

ALPHA  -  GAMMA) * (U/2 . ) 

BETA  -  THETA/2 .)* (DPRIME/2 . ) 

* (XLENl  +  XLEN2) 


Q********** ******  NOW  CORNER  IS  COMPUTED  ******************************* 


IF  (U  .LE.  WO)  THEN 

AREAl  =■  .  5*PI/3  .  *  (  (PD  +  2  .  *WO  -  U)/2.)**2 
AREA2  =  . 5*PI/3 .* (DPRIME/2 .) **2 
ECZ  =  AREAl  -  AREA2 

XLENl  =  PI/3.* ((PD  +  2.*WO  -  U) /2 . ) 

XLEN2  =  PI/3 .*  (DPRIME/2 . ) 

SCZ  =  (GL  -  U)* (XLENl  +  XLEN2 ) 

ELSE 

ECZ  =0.0 
SCZ  -  0.0 
END  IF 


Q**  ***********  A- I'r  *  *  FINALLY  THE  INNER  SLIVERS  ****★*****•******■*•■**•**•*■★■**  *  * 

(2k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  it  k  k  k  k  k  k  k  k  k  -k  k  k  k  -k  -k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  -k  .'t  k  k  k  k  k  if  k  k  k  k  k  k  k  -k  k  k 


SQRT3  =  SQRT 
TOPP  “  (WI 
IF  (U  .LE. 
AREAl  = 
AREA2  = 
AREA3  = 
EIZ  " 
XLENl 
SIZ  = 

ELSE 


(3. ) /2. 

+  PD) /SQRT3  -  PD 
TOPP)  THEN 
=  . 5*SQKT3*  (WI  +  PD) **2 
=  .5*(W1  +  PD) * (DPRIME/2 .) *SIN 
=  .5* (PI/3,  -  THETA) * (DPRIME/2 , 
AREAl  -  3.*(AREA2  +  AREA3) 

=  (PI/3  -  THETA) * (DPRIME/2 . ) 

3 . * (GL  -  U) *XLEN1 


(THETA/2 . ) 
)  **2 


EIZ  =  0.0 
SIZ  =0.0 


END  IF 


Q^kkkkkkkkkkkkkkkkkkkkkk'kkkkkkkkir-Akkkkkkkkkkkkkk''{kkkkkkkkkkn'tKk  t^Xkkkkkkkk-k 

C****************  SURFACE  AREA  AND  UNBURNED . VOLUME  COMPUTED  ************ 

(;;***  ***Tr**************************A*  *-***•  ************************  ********* 


E  NIZ*EIZ  +  NOZ*EOZ  +  NCZ*ECZ 

SFAREA  =  2 . *E  +  NIZ*SIZ  +  NOZ*SOZ  ^  NCZ*SCZ 

VULUNB  -  E* 'GL  -  U) 

RETURN 
END  IF 

Q^k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  -k  -k 

(2k  k  k  k  k  -k  -k  k  -h  k  k  k  A  k  k  k  k  k  k  k  k  k  p^QliLED  BALL  GRAIN  kkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

Qkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk'kkickkkkirkkkk 

IF  (TCODE  .EQ.  14)  THEN 
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DD  “  A3  (8)  -  A3  (7) 

TT  ■=  A3  (7)  -  U 

SFAREA  =  PI/2, *(DD**2  +  2.*TT**2  +  PI*DD*TT) 
VOLUNB  =  TT*PI/2 . *  (  (DD**2) /2 .  +  (TT**2)/3.  + 
RETURN 


(PI*DD*TT) /4 . 


END  IF 


-K-k-k-kk-k-k-kic-k-ki 


Q-k  k  k  *  -k  k  k  k  -k  k  -k  k  k  k  k  k  k  k  -k  k  k  -k  CIGARETTE  GRAIN  ************************************ 

0-AA-*i>r*-******A*  A  ArAA-********  **********************  *'*****it***-A:*)»r  A  *★★★*■**•**** 

IF  (ICODE  .EQ.  13)  THEN 
SFAREA  =  PI*D*D/4 
VOLUNB  -  (GL  -  R)*SFAREA 
RETURN 
END  IF 

(^k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

C************** ***  OTHER  GRAINS  ARE  HANDLED  ************************ 

GO  TO  (  1000,  1010,  1020,  1030,  1040,  1050,  1060,  1070, 

1  1080,  1090,  1100,  1110),  ICODE 

(^k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  rr  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

Qk k k k k k k k h k k k k k k k k k k  CODE  ll  7“ PERF  GRAIN  k k k k k k k k k k k k k k k k k k k k k k k k k k k k 

(2'*'  k  k  k  k  k  K  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

C***  This  part  calculat.es  the  conditions  before  the  grain  burns. 

Qk  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  \  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

1050  CONTINUE 

D  =  3 . 0*PD  r  2 . 0* (WI  +  WO) 

EO  =  PI*(D**2  -  7 .0*PD**2) /4 .0 
SO  =  PI*(D  +  7.0*PD)*GL  +  2.0*E0 
VO  =  E0*GL 
WW  -  WI  +  PD 
DO  1120  K  -  1,  3 
S7  (K)  =■  WW 
1120  CONTINUE 

WEBC  -  AMINl  (WO,  WI,  GL) 

QkkkkkkkkkkkkkkitkkkkkkkkkkitkkkkAkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

C***  This  part  does  ttie  calculations  for  the  burning  grain. 

Qkkkkkkkkkkkkkirkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

GRL  =  AMAXl  (GL  -  U,  0.0) 

OD  =  D  -  U 
PHFD  ==  PD  +  U 

IF  (U  .GT,  WEBC)  GO  TO  1130 

E  =  PI*(OD**2  -  7 .0*PRFD**2) /4.0 

SFAREA  -  PI*(OD  +  7 . 0 *PRFD) *GRL  +  2.0*E 

FRCSFA  =  SFAREA/SO 

VOLUNB  =  E*GRL 

VOLBRN  =  VO  -  VOLUNB 

FRCBRN  -  VOLBRN/VQ 

RETURN 

g********* ******* ************A***** ************************  ************** 

C***  This  part  does  the  calculations  for  when  the  grain  slivers.  ****** 
g***********************  *****************************  ***************.***•*. 
1130  CONTINUE 

CALL  GENIS  (S7,  PHFD,  GRL,  SFl,  GVl) 

CALL  GENOS  (S7,  PRFD,  GRL,  0.5*OD,  SF2,  GV2 ) 


•kkkkkkkkkkkkkkkkkkkkkkkkkk 


:  k  k  k  k  k  k  kkkkkkkk 
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SFAREA 

FRCSFA 

VOLUNB 

VOLBRN 

FRCBRN 

RETURN 


=  6.0* (SFl  +  SF2) 
=  SFAREA/SO 
=  6.0* (GVl  +  GV2) 
=  VO  -  VOLUNB 
-■=  VOLBRN/VO 


ii:  if  i(  -k  if  "k  if  "k  ir  k  -k  i!  ic  -k  -k  -k  -k  "k  "k  "k  -k  "k  it  k  k  ie  -k  -k  ic  ic  i(  ic  ic  ie  ic  ic  ic  ic  'k  ic  ic  “k  'k  ir  ir  "k  -k  "k  -k  'k  ic  if  if  -k  "k  -k  ic  -k  -k  -k  -k  -k 

Q'k'k -k-k -k -k -k-k -k -k -k -k -k -k -k -k  *  ^  it  CODE  2  I  1~PERF  GRAIN  kkkkkkkkkkkkkkkkkkkkkkkkkkkk 

^kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk-kk'kk'kkkkkkkkkkkkkkkkkk'kkkkkkkkkkkkkkkkk 

Q***  This  part  calculates  the  conditions  before  the  grain  burns. 

(^kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk'kkkkk 


1030  CONTINUE 

D  -  PD  2.0*WI 
EO  =  PI*(D**2  -  PD**2)/4.0 
SO  =  PI*(D  +  PD)*GL  +  2,0*E0 
VO  =  E0*GL 

WEBC  =  AMINI  (GL,  WI) 

Q^kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk'iek-k'kk'kkk’kickk'k'kkk'kk'k'kkk'kk 

C*******  This  part  does  the  calculations  for  the  burning  grain.  ******* 


IF  (U  .GE.  WEBC)  TtlEN 
GRL  =0.0 
E  ■  0.0 
GO  TO  IHO 
END  IF 

GRL  =  GL  -  U 
OD  =  D  -  U 
PRFD  =  PD  a  U 

E  =  PI*(OD**2  -  PRFD**2)/4.0 
1140  CONTINUE 

SFAREA  =  PI*(OD  +  PRFD) *GRL  +  2.0*E 

FRCSFA  =  SFAREA/SO 

VOLUNB  =  E*GRL 

VOLBRN  =  VO  -  VOLUNB 

FRCBRN  =  VOLBRN/VO 

RETURN 

(^irkkkkkkkkkkkkkkkkkkkkk  kkkkkkkkkkkkkkkkkkkkkkkkkkkk  kkkkkkkkkkkkkkkkkkk-^  k 

C*******. ■.*********•**  CODE  J;  CORD  GRAIN  ****************************** 

(^■kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk'kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

C***  This  part  calculates  the  conditions  before  the  grain  burns.  ****** 

Qk  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  i:  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  -k  k 

1010  CONTINUE 

SO  =  GL*PI*D  +  PI*D**2/2.0 
VO  =  GL*PT*D**2/4 . 0 

C***  This  part  does  the  calculations  for  the  burning  brain. 

GRL  -  AMAXl  (GL  -  U,  0.0) 

OD  -  A1-1AX1  (D  -  U,  0.0) 

E  =  PI*OD**2/4 .0 

SFAREA  =  PI*OD*GRL  +  2.0*E 

FRCSFA  =  SFAREA/SO 

VOLUNB  =  GRL*PI*OD**2/4  ,0 

VOLBRN  =  VC  -  VOLUNB 

FRCBRN  =  VOLBRN/VO 

RETURN 


2« 


(^**************  A*  **********<r******-***-Jf*  *****-*******  A- ****■*>-  A' A-  A- 

C*****A*********A:***  CODE  4:  RECTANGULAR  STRIP  GRAIN  **************** 

Q***it*AArA*Ar*AA**AAr*A’7tA’Tlr‘A‘A-A‘A*Ar****Ar***’Ar'Jlr’A*-‘A-A****Ar****-A:*-A-A**A-*A**A'-AA:A'AA*-*'* 

.1020  CONTINUE 

SO  =  2.0*(GL*D  +  D*WI  +  GL*WI) 

VO  =  GL*D*WI 

GRL  =  AMAXl  (GL  -  U,  0.0) 

DS  -  AMAXl  (D  -  U,  0.0) 

WIS  =  AMAXl  (WI  -  U,  0.0) 

SFAREA  -  2.0*(GRL*DS  +  DS*WIS  +  WIS*GRL) 

FRCSFA  -  SFAREA/SO 
VOLUME  =  GRL*DS*WIS 
VOLBRN  =  VO  -  VOLUME 
FRCBRN  =  VOLBRN/VO 
RETURN 


C* 


■k  :k  -A  ■kk-kk-kkkkkkkkifAfirkk'k  • 


*  CODE  5:  SPHERICAL  GRAIN 


r-Mkkkirkkk'M'ki 


*****ArAAAAr*AAr* 
■  kkk  k  k 
r-kkkkk’k'k'k'k'kAkk 


1000  CONTINUE 

50  -  PI*D**2 

VO  =  PI*D**3/6.0 

OD  =  AMAXl  (D  -  U,  0.0) 

SFAREA  =  PI*OD**2 
FRCSFA  -  SFAREA/SO 
VOLUME  ■-  Pl*OD**3/6 . 0 
VOLBRN  ■-=  VO  -  VOLUME 
FRCBRN  -  VOLBRN/VO 
RETURN 

Qk  k  A  A  k  k  k  k  A  k  k  k  k  k  k  A  k  A  k  k  A  k  k  k  -k  A  k  k  k  A  k  k  k  -k  k  k  k  k  k  k  k  k  k  k  -k  k  A  k  A  A  A  A  A  A  k  A  k  k  A  A  k  k  A  k  k  A  k  k  k  A  k 

QkkkAkkkkkkkkkAkkkk  CODE  6;  SLOTTED—TUBE  GRAIN  *•  *  *  *  *  *  5*r  A- A- tA  Ar  A  *  A -A  A- A  A  ★  *  Ar  A 

QkkkkkkkkkkAkkkkkkkkkkkkkkkkkkkkkkAkkAkAkkkkkkAkkkAkkkAkkkkAAkAkkAkkkkkk 

C***  This  part  doas  the  calculations  before  the  grain  burns.  ********** 

qaaaaaaaaaaaaaaaaaaaaaaa.vaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa-aaaaaaaa 

1040  CONTINUE 

SLOT  -  0.5*WM 
so  =  0 . 5*D 

51  =  0 . 5*PD 

THETA  =•  ASIN  (SLOT/SO) 

ALPHA  =  ASIN  (SLOT/SI) 

EO  =  (PI  -  ALPHA) * (SO**2  -  SI**2)  +  (SO  -  S I ) * *  2 * ALPHA 

SO  •*  2.0*  ((PI  -  ALPHA)  *SI  +  (PI  -  THETA)  *SO  +  (SO*COS  (THETA) 

1  -  SI*COS  (ALPHA) )) *GL  +  2.0*E0 

VO  =  GL*E0 
WI  F=  SO  -  SI 
WEBC  =  AMINl  (GL,  WI) 

(^AkkAkkkkk  A  kkkkkkkkk  A  kAAkkkkkk  A  AkAkkAkkAkkAkk  A  kkkkkkkkkk  A  A  kkkkkkkAkkk  A  kk 

C*A*  This  part  does  the  calculations  for  the  burning  grain,  a********** 

(2^AkA><kkkAkAkkAAAkkAAkkkkkkkAkAkkkAkAAAkkkkkAkAAkkAAAkAAkkAkAkkAAkAKkAAAk 


IF  (U  .GF,  ,  WEBC)  THEN 
GRL  =0.0 
E  -  0.0 
GO  TO  1150 
END  IF 

SLOT  -  0 . 5* (WM  +  U) 


so  =  0 . 5*  (D  -  U) 

SI  =  0.5*  (PD  +  U) 

GRL  =»  GL  -  U 

THETA  -  ASIN  (SLOT/SO) 

ALPHA  =  ASIN  (SLOT/SI) 

E  =  (PI  -  ALPHA) * (SO**2  -  SI**2)  +  (SO  -  SI)**2*ALPHA 
1150  CONTINUE 

SFAREA  -  2.0* ((PI  -  ALPHA) *ST  +  (PI  -  THETA) *SO  +  (SO*COS  (THETA) 
1  -  SI*COS  (ALPHA) )) *GRL  +  2.0*E 

FRCSFA  -  SFAREA/SO 
VOLUNB  =  E*GRL 
VOLBRN  =  VO  -  VOLUNB 
FRCBRN  =  VOLBRN/VO 
RETURN 


(2*'A’*****  +  ^Tlr'*'**A  kii  kk-kk'kic-k^kit-kiiic'kk'k'k-k-k-k-k-kic-kiKie-kicic-kitiK’k-k'k-kir-k-ic-kick'^'k'k'ic-k'k-kiririKir-k*: 

Q** ***********  A ** ***  CODE  7:  ROUND-HEX  37-PERF  GRAIN  ***************** 

(^kkkkkkkkkkk'kkkkkkkiKk-k'kkkkkk-k'k'k-k'k'kir-kic-k-Ak-kk-k'kic-k-k-k-k-ki^k-k-k-k-kif-kk-kkk'k-k'k'kk-k-kk-k-k 


1090  CONTINUE 

50  =  18  . 

51  =  54  . 
NPERF  =  37 


D  =  7.0*PD  +  6.0*WI  4  2.0*WO 
GO  TO  1160 


C** **************************************************************  ******* 

C** . ^ ************ ***  CODE  8;  ROUND-HEX  19-PERF  GRAIN  ***************** 

Q****'*-***^:^*-*  *******-**/rA'Ar*yt:ArA-***-*  A:  **•*■* -A- ir'Ar**'*-*********-A'-ilryr*'Ar7S:**'Ar**  *•★*■*★★★*: 


1080 


CONTINUE 

50  =  12.0 

51  »  24.0 


NPERF  =19 

D  -  5.0*PD  +  4.0*W1  +  2.0*WO 

Qkkkkkie'k-k'k’kif'kk-kk-kkk’k’kirieirkkk'kkkkk'k'kic-kit-kfcii'kkic-k-k’k'k'k-kk'k'kkirkir-kk-kkk'. 

C******************  CALCULATIONS  FOR  CODES  7,8,10 

Qk  k  k  k  k  k  -k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  A  k  k  k  k  k  : 

1160  CONTINUE 

WW  =  VJI  +  PD 
WW2  =  WW**2 


;  A'  *  ★  ★  -A:  ★  A:  Ar  A  k  k  k  k  k  k  k  k  k  k  i: 


PRFD  ■  ■  PD  (•  U 

PRFD2  =  PRFD* *2 

GRL  =  AMAXl  (GL  -  U,  0.0) 

E  -=  0 . 0 

THETA  =  2.0*ACOS  (AMINl  (WW/PRFD,  1.0)) 

ALPHA  =  ACOS  (AMINl  ((2.0*WO  4  PD  -  U)/PRFD,  1.0)) 

IF  (U  .LT.  WO)  E  =  0 . 25*PI*  (  (2 .0*WO  +  PD  -  U)**2  -  PRFD2) 

IF  (THETA  ,GE.  PI3)  GO  TO  1170 

E  =  E  4  SI*0 . 25* (WW2*RT  -  1 . 5*PRFD2 * (SIN  (THETA)  4  PI3  -  THETA)) 
1170  CONTINUE 

IF  (ALPHA  .GE.  0.5* (PI  -  THETA))  GO  TO  1180 
E  =  E  4  SO*0 . 125* (2 . 0* (2 . 0*WO  4  PD  -  U)*(2.0* 

1  WW  -  PRFD*SIN  (ALPHA) ) 

2  -  PRFD2*(SIN  (THETA)  ^  PI  -  2.0*ALPHA  -  THETA)) 

1180  CONTINUE 

IF  (2.0* WO  4  PD  .LT.  Wl)  THEN 

WRITE  (*,  *)  ' *FORMT*  BAD  HEX  PROP' 
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PAUSt: 

GO  TO  1190 
END  IF 

VOLUNB  =  E*GRL 

C*******************  xEST  TO  SEE  IF  GRAIN  CONSUMED  ******************** 
IF  (VOLUNB  .LE.  0.0)  THEN 
SFAREA  ---  0.0 
VOLUNB  -0.0 
GO  TO  1190 
END  IP 

A  A-A-AA-A:  THE  S  UREACE  AREA 

PH  -  D/2.  -  wo  -  PD/2. 

IF  (U  .EQ.  0.0)  THEN 

SFAREA  =  2  .  *•£  i-  GRL*PH*6.  +  NPERF*PI*PD* 

I  GRL  +■  PI*GRL*(2*WO  +  PD) 

GO  TO  1190 
END  IF 

QAAAAAAAAAAtAAA-AAAAAAAAAAAAA  j\JQ  S  El  VER  TNG  YET***’^***^’''*’''*  AAtA^'A:*  AAAA*AAA*A 
IF  ((WO  .GT.  U)  .AND.  (WI  . GT .  U) )  THEN 

SFAREA  -  2.*E  +  NPERF* (PD  +  U)*GRL*PI  +  6.*PH* 

1  GRL  +  PI*(2.*WO  +  PD  -  U) *GRL 

GO  TO  1 1 90 
END  IF 

(^ ******************* *  SLIVERING  ********************************** 

SFAREA  2.*E 

(^****** *****************  fxR,ST  the  inner  SLIVERS  ********************* 

IF  (THETA  .GE.  PI3)  THEN 
GO  TO  1200 
ELSE 

SFAREA  “  SFAREA  *  1 . 5 *PRFD*GRL* ( PI 3  -  THETA) *SI 
END  IF 

1200  CONTINUE 

e***************  NOW  THE  OUTER  SLIVERS  &  CORNERS  ********************** 
0* ****************** *  CORNERS  NOT  CONSUMED  **************************** 
IF  (WO  .GT.  U)  THEN 

SFAREA  -  SFAREA  PI*  (PD  H-  2 .  *WO  -  U  PRFD)  *GRL 
END  IF 

0************...***  **  ***NOW  OUTER  SLIVERS********  **********  **  ***  *  ****** 
IF  (ALPHA  .LT.  .5*  (PI  -  THETA))  THEN 

SFAREA  -  SFAREA  h  (WW  -  PRFU*SIN  (ALPHA) ) *GRL*SO  + 

1  PRFD*GRL*  (PI/2  .  -  ALPHA  -  TllETA/2  .  )  *SO 

END  IF 

1190  CONTINUE 
RETURN 

^-kk^f^-k-k'k'k-k'^-k'Ak-khk'^ii^'h'k'h'k'k'kkiticicirfeicAicAk'.  ff-k-k-kicic-k-kiv-k'k'k-k’kA-k'kirir-k-k'k'K'M'kAM'kk^'k'k-kic 

Q* -k -k  ^ -k -k -k -k -k -k -k -h -K  k  k  k  -k -k -k  CODE  9:  19  — PERF  GRAIN  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

(j-kkkkkkkkkkkkk**,kk-h\ktr-k;c-icirA^kkkkk*:kkkkkkkkiikkk*kkkk-kfzkkkkkkkkkkkkkkt.  kkkk 

C***  This  part  calculates  the  conditions  before  tlie  grain  burns  ******** 

Q*******'******A***:A:***:Ar*-*-*******-******A*A:*  ******************************** 

1070  CONTINUE 

D  -  b.0*PD  (-  2,0*  (Wl  +  WM  +  WO) 

EO  =  PI*(D**2  -  19 . 0*PD**2) /4 .0 
SO  ■=  PI*(D  *  19.0*PD)*GL  *  2.0*E0 
VO  -  E0*GL 


S19 

(1, 

1) 

= 

WI  +  PD 

G19 

(2, 

1) 

S19  (1,  1) 

S19 

(3, 

1 ) 

S19  (1,  1) 

Sx9 

(1, 

2) 

= 

0.5*SQRT  (3.0* (WM 

S19 

(2, 

2) 

= 

S19  (1,  2) 

S19 

(3, 

2) 

=: 

S19  (1,  1) 

S19 

(1, 

3) 

sr 

PD  I  0 . 5* (WI  +  WM) 

S19 

(2, 

3) 

S19  (1,  2) 

S19 

(3, 

3) 

= 

WM  +  PD 

S19 

(1, 

4) 

SI  9  (1,  3) 

S19 

(2, 

4) 

= 

2.0*S19  (1,  3) 

S19 

(3, 

4) 

= 

SI  9  (1,  3) *RT 

WEBC 

= 

AMINl 

(WO,  WM,  WI,  S19  ( 

+  PD) **2  +  (WI  +  PD) **2) 


(1,  3)  -  PD,  S19  (1,  2)  -  PD,  GL) 
C***  This  part  does  the  calculations  for  the  burning  grain.  ************ 

Q**:****-^-*:*-*-*-***'*-*'***-*-****'*-***^*  *•***•*■*#****••**★★*:*  yi-**7kr**  >*•*★•*•  -A- ****■*•**★★★** 

GRL  =  AMAXl  (GL  -  U,  0.0) 


OD  =  D  -  U 
PRFD  =  PD  +  U 

IF  (U  .GE.  WEBC)  GO  TO  1210 
E  =  0 .25*PI* (OD**2  -  19 . 0*PRFD**2) 

SFAREA  -  PI*(OD  +  1 9 . 0*PRFD ) *GRL  +  2.0*E 

FRCSFA  =  SFAREA/ SO 

VOLUNB  =  E*GRL 

VOLBRN  -  VO  ~  VOLUNB 

FRCBRN  ’  VOLBRN /VO 

RETURN 

1210  CONTINUE 


•SUMS A  -  0.0 
SUMGV  =0.0 
DO  1220  K  =  1,  2 

CALL  GENIS  (S19  (1,  K) ,  PRFD,  GRL,  SA,  GV) 
SUMSA  =  SUMSA  +  6.0*SA 
SUMGV  =  SUMGV  +  O.O^GV 
1220  CONTINUE 

CALL  GENIS  (S19  (1,  3),  PRFD,  GRL,  SA,  GV) 

SUMSA  =  SUMSA  +  12.0*SA 
SUMGV  =  SUMGV  +  12.0*GV 

CALL  GENOS  (S19  (1,  4),  PRFD,  GRL,  .5*00,  SA,  GV) 

SUMSA  =  SUMSA  +  12.0*SA 

SUMGV  =  SUMGV  +  12.0*GV 

SFAREA  SUMSA 

FRCSFA  -  SFAKEA/SO 

VOLUNB  SUMGV 

VOLBRN  =  VO  -  VOLUNB 

FRCBRN  =  VOLBRN/VO 

RETURN 


k'kkA'k^AA'k  A'  ^'kk-k'k-k-icieic-k-ktf'frlc'k^k'k'A'k'kir'k'k-k-k'kkk'kkk  if  Mifk'k'k'fc'kkir'k'kif 

C*******************  20;  ROUND-HEX  7-PERF  GRAIN  ***************** 

(''■^■kA-k^ir-kk^k^k-k-k-kk^ic-kk^kifk'kk'k'kk-kifit-kfriciiir-k'kkkkirick'kkAk'k'kifk'klKirk-kkici^ickkkk-kKic'k'fr 


1060  CONTINUE 
SO  =  6.0 
Si:  6.0 
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Vilbat  RknUniA.  WariUikliM  A>Tb . 


NPERF  -  7 
GO  TO  1160 

^irieici^iriciriririfie'kic'kkif'kic'k'k'kicicir'kir-k-k'k'k'k'k'k-k'kic-k-k’kit-kicir’k'ic'kieic'kifitic-kitic-k'k-k'k-ic-k-k-ic^-k-k-k-^-k-k-ie 

C*******************  CODE  11:  CORD  WITH  INHIBITED  ENDS  ***************** 

^•k'kic-k-Air'k-kif-k-A  A  ic^'M-ki^-k-kir'k'k'k-kieiric-k-k-kif-kif-kick-kir-kie'k'k-kir ‘Jr  ic'k'kif-k-kififirii  A  ■kk’kic'kic'k'k'k-kic'AA'kk 

C  This  routine  will  only  calculate  the  surface  area  of  the  lateral  surfa 

Cces . 

C  It  will  not  calculate  the  surface  area  of  the  inhibited  ends.  ******** 

Q*  *  *  ***************  k  *  *  *  *  -k  ******  -k  k  ************************************** 

1100  CONTINUE 

SO  =  GL*FI*D 

VO  -  GL*PI*D**2/4  .  0 

C  ****  This  part  does  the  calculations  for  the  burning  grain 
OD  -■=  AMAXl  (D  -  U,  0.0) 

SFAREA  ■  PI*OD*GL 
FRCSFA  -  SFAREA/SO 
VOLUNB  =  GL*PI*OD**2/4 . 0 
VOLBRN  =■  VO  -  VOLUNB 
FRCBRN  =  VOLBRN/VO 
RETURN 

Q*********v  ★*★*!★**★★★★**★**********★****★★*★★★★*■**★★  •*★'***’*★★*★★★★*★★★★★★ 

C** ********  CODE  12:  RECTANGULAR  STRIP  GRAIN  WITH  INHIBITED  SIDES  ****** 

^•k  -k  ie  ir  -k  if  "k  -k  ir  -k  -k  "k  -k  if  -k  ‘k  ir  -k  -k  k  -k  ■k  ^  -k  -k  ir  -k  k  -k  ir  k  -k  -ic  "k  "k  if  ic  ic  1C  -k  k  -k  ir  -k  ir  i(  ^  -A  -k  -k  ic  -k  -A  "ic  if  if  ^  it  -k  -k  -k  ir  ^  "k  if  "k  it  "k  i< 

C  This  routine  will  only  calculate  the  surface  area  of  the  two  burning  s 

Cides . 

C  It  will  not  calculate  the  surface  area  of  the  inhibited  sides. 

Q**^t#**’*Tlr*  +  **»**>^*A’*’Tlr***-*r****'*-*'*r-*************'**-jir**-*******'**-*-^*****'Ar'***** 

1110  CONTINUE 

SO  =  2 . 0*GL*WI 
VO  =  GL*-D*w; 

VOLMAO  =  (D  -  WI)*GL*WO 
VOLMBO  =  WI*GL*WO 

C  ***  This  part  does  the  calculations  for  the  burning  grain 
DS  =*  AMAXl  (D  -  U,  0.0) 

SFAREA  =  2.0*GL*WI 
FRCSFA  -  SFAREA/SO 
VOLUNB  =  GL*DS*WI 
VOLBRN  =  VO  -  VOLUNB 
FRCBPN  =  VOLBRN/VO 
IF  (D.,  .GE.  WI)  THEN 
VOLABR  =  VOLBRN 
VOLBBR  =  0 . 

ELSE 

VOLABR  =  VOLMAO 
VOLBBR  =  VOLBRN  -  VOLABR 
END  IF 
RET’.  IN 
END 

0*****A****'*  ■*■  ************:  ****-***********A:Ac***  ****************  ★*■*■*******■*■ 

C 

C  SUBROUTINE  *GENIS* :  calculate  surface  area  and  volume  for  a 

C  general  inner  sliver  of  a  burning  grain 
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oonnn  nnno  oooo  onnn  non  nonno  on 


with  length  =  GRL  &  perforation  dia.  =  PRFD. 

SUBROUTINE  GENIS  (S,  PRFD,  GRL,  SURF,  VOL) 

DIMENSION  S  (3),  A  (4) 

DATA  PI2/1 . 5707963/ 


**********  :  store  angles  A1,A2,A3  and  area  of  triangle 

with  sides  S (1) , S (2) , S (3)  into  A(1)...A(4) 

A  (1)  -  ACOS  ((S  (2)**2  +  s  (3)**2  -  S  ( 1) * *2 ) / ( 2 . 0 *S  (2)*S  (3))) 

A  (2)  =  ACOS  ({S  {1)**2  +  S  (3)**2  -  S  (2) **2) / (2 . 0*S  (1)*S  (3))) 

A  (3)  =  ACOS  ((S  (1)**2  +  S  (2)**2  -  S  (3 ) **2 j / (2 . 0 *S  (1) *S  (2))) 

A  (4)  =  0.5*S  (1)*S  (3)*SIN  (A  (2)) 

...check  for  error  condition:  find  if  triangle  acceptable... 

J  =  0 

DO  1000  I  =  1,  3 

IF  (A  (I)  .LT.  0.5*PI2)  J  =  J  +  1 
1000  CONTINUE 

IF  (J  .GT.  1)  STOP  '  GENIS  ERROR' 


succeeding  passes  until  burnout:  find  auxiliary  angles 

TAU12  =  ACOS  (AMINl  (1.0,  S  (3) /PRFD)) 

TAU13  -=  ACOS  (AMINl  (1.0,  S  (2) /PRFD)  ) 

TAU23  =  ACOS  (AMINl  (1.0,  S  (1) /PRFD) ) 

...and  branch  to  25  if  sliver  fails  burnout  criteria.  If  not 

then  sliver  is  burned  and  go  to  30. 

IF  (TAU12  +  TAU13  +  TAU23  . LT .  PI2  .AND.  GRL  . GT .  0.0)  THEN 
GO  TO  1010 
ELSE 

GO  TO  1020 
END  IF 


sliver  not  burned  out:  determine  end  area,  volume  and  surface  area 
1010  CONTINUE 

E  =  A  (4)  -  0.25*PRFD*(S  {1)*SIN  (TAU23)  +  S  (2)*SIN  {TAU13) 

1  +-  S  (3)*SIN  (TAU12)  e  PRFD''  (PI2  -  TAU12  -  TAU13  -  TAU23)) 

VOL  =  E*GRL 

SURF  =  2.0*E  +  GRL*PRFD* {PI2  -  TAU12  -  TAU13  -  TAU23) 

. . . and  RETURN 

RETURN 

C 
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n  n  c  n  n  n  o  o  n  c  o  o  n  o  o  n  n  n  o  n  n  n  n  n  n  no 


sliver  is  burned  out:  return  with  zero  volume  and  surface  area. 


1020  CONTINUE 
VOL  =0,0 
SURF  =  0.0 
RETURN 
END 

*Tk-*-VV+*-<r-’^+-*-*  +  k-+i»A-*A:*A*****A*A***r*******-****?ir*****A**Tlr-*********'******** 

SUBROUTINE  "GENOS"  :  Calculates  surface  area  and  volume  for  a 

general  outer  sliver  of  a  burning  grain 
with  length  =  GRL,  radius  =  RAD,  and 
perforation  diameter  =  PRFD 


SUBROUTINE  GENOS  (S,  PRFD,  GRL,  RAD,  SURF,  VOL) 
DIMENSION  S  (3) ,  A  (4) 


•k'Aie'A'Ai.'-kir-AJr-kA 


Store  angles  A1,A2,A3  and  area  of  triangle 
with  sides  S  (1) , S (2) , S  (3)  into  A(l)  ...A(4) 


A  (1)  =  ACOS  ((S  (2)**2  +  S  (3)**2  - 

A  (2)  =  ACOS  ((S  (1)*^2  +  S  (3)**2  - 

A  (3)  =  ACOS  {(S  (1)**2  +  S  (2)**2  - 

A  (4)  =  0,5*3  (1)*S  (3)*SIN  (A  (2)) 


S  (1) **2) / (2 . 0*S  (2) *S  (3) ) ) 
S  (2) **2) / (2 . 0*S  (1)*S  (3))) 
S  (3) **2) / (2 .0*S  {1)*S  (2))) 


succeeding  passes  until  burnout;  determine  auxiliary  angles 

TAUl  =  ACOS  (AMINl  (1.,  (S  (2)**2  +  RAU**2  -  0.25* 

1  PRFD**2) / (2 . *S  (2)*RAD))) 

TAU2  =  ACOS  (AMINl  (1.,  (S  (3)**2  +  RAD**2  -  0.25* 

1  PRFD**2) / (2 . *S  (3)*KAU))) 

TAU3  =  ACO.S  (AMAXl  {  -  1.0,  (S  (2)**2  -  RAD**2  +  0.25* 

1  PKFD**2)/(S  (2)*PRFD))) 

TAU4  =  ACOS  (AMAXl  (  -  1.0,  (S  (3)**2  -  RAD**2  +  0.25* 

1  PRFD**2)/(S  (3)*PR.FD))) 

SIG  =  ACOS  (AMINl  (1.0,  S  (D/PRFD)) 


.  .  .then  check  error  conditions.  . . 

IF  (TAU3  .LT.  A  (3)  .OR.  TAU4  . LT .  A  (2))  STOP  '  *GENOS*  ERROR' 


IF  ok,  check  if  sliver  burned  out.  If  not  burned  out  go  to  25. 
If  burned  out  go  to  30. 

IF  (TAUl  I  TAU2  .  LT .  A  (1)  .AND.  GDI.  .  GT .  0.0)  THEN 
GO  TO  1000 
ELSE 

GO  TO  1010 
END  IF 


sliver  not  burned  out;  determine  end  area,  volume  and  surface  area. 
C 


35 


CJ  CJ  o  u  u  u 


1000  CONTINUE 

E  =  0.5*RAD*(S  (2) ‘STM  (TAUl)  +  RAD* (A  (1)  -  TAUl  -  TAU2) 

1  +  S  (3)*SIN  {TAU2))  -  A  (4)  -  0.25*PRFD*{S  (1)*SIN  (SIG) 

2  +  0 . 5*PRFD* (TAU3  +  TAU4  -  2.0*SIG  -  A  (2)  -  A  (3))) 

C 

VOL  =  E*GRL 
C 

SURE  -  2.0*E  +  GRL*(RAn*{A  (1)  -  TAUl  -  TAU2)  +  0 . 5 *PRFD* (TAU3 
1  +  TAU4  -  2.0*SIG  -  A  (2)  -  A  (3))) 


.  . . and  RETURN . 
RETURN 


sliver  is  burned  out:  return  with  zero  volume  and  surface  area. 

1010  CONTINUE 
VOL  =0.0 
SURF  =0.0 
RETURN 
END 
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U.S.  Army  Armament  Research,  Development, 
and  Engineering  Center 
AITN:  SMCAR-QAH-T,  Mr.  John  Domcn 
Building  62  North 
Pieatinny  Arsenal,  NJ  07806-5(XX) 

1  Commander 

U.S.  Army  Armament  Research,  Development, 
and  Engineering  Center 
ATTN:  SMCAR-AEE-ER, 

Mr.  Ken  Klingiiman 
Building  1501 

Pieatinny  Arsenal,  NJ  07806-5(XX) 

1  Commander 

U.S.  Army  Annainent  Research,  Development, 
and  Engineering  Center 
ATTN:  SMCAR-ASI,  Dr.  Jerry  Rubin 
Pieatinny  Anscnal,  NJ  078()6-5(XX) 

1  Director 

U.S.  Army  BMD 
Advanced  Technology  Center 
P.  O.  Box  1500 
Huntsville,  AL  35807 

1  Chairman 

DOD  Explosives  Safety  Board 
Room  856-C 
Hoffman  Bldg.  1 
2461  Eisenhower  Ave. 

Alexandria,  VA  22331-0600 

1  Department  of  lire  Army 

Office  of  the  Product  Manager 
155mm  Howitzer,  M109A6,  Paladin 
ATTN:  SFAE-AR-HIP-IF,  Mr.  R.  De  Kleinc 
Pieatinny  Arsenal,  NJ  07806- 5(XX) 

1  Commander,  USACECOM 

R&D  Technical  Library 
A'rTN:  ASQNC-ELC-IS-L-R,  Mycr  Center 
Fort  Monmouth,  NJ  07703-5301 


3  PEO-Amiamcnts 
Project  Manager 
Tank  Main  Armament  Systcm.s 
ATTN:  AMCPM-TMA,  K,  Russell 
AMCPM-TMA-105 
AMCPM-TMA-i20 
Pieatinny  Arsenal,  NJ  07806-.5(XX) 

1  Commander 

Production  Base  Modernization  Agency 
U.S.  Army  Amiaincnt  Rcscareh, 
Development,  and  Engineering  Center 
ATTN;  AMSMC-PBM-E,  L.  Laibson 
Pieatinny  Arsenal,  NJ  07806-5000 

5  Director 

Benet  Laboratories 

U.S.  Army  Walervlict  Arsenal 

ATTN:  SARWV-RD, 

G.  Caralano 
R.  Thierry 
R.  Hasoenbein 
P.  Vulis 
P.  Alto 

Waicrvliei,  NY  12189 

3  Commander 

U.S.  Army  AMCCOM 
A'lTN;  AMSMC-IRC,  G.  Cowan 
.SMCAR-ESM(R), 

W .  Fortune 
R.  Zasmow 

Rock  Island,  IL  61299-7300 

1  Commandant 

U.S.  Army  Aviation  Sehtx)! 

ATTN:  Aviation  Agency 
Fort  Rucker,  AL  36360 

1  Director 

HQ,  TRAC  RPD 

ATl’N;  ATCD-MA,  MAJ  Williams 
Fort  Monrtre,  VA  23651-5143 

1  Headquarters 

U.S.  Army  Materiel  Command 
An'N;  AMCICP-AD,  Michael  F.  Fisette 
5(X)1  Eisenhower  Ave. 

Alexandria,  VA  22333-(XX)l 
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4  Commander 

U.S,  Army  Armament  Research, 

Development,  and  Engineering  Center 
ATTN;  SMCAR-CCD,  D.  Spring 
SMCAR-CCS 

SMCAR-CCH-T,  L.  Rosendorf 
SMCAR-CCH-V,  E.  Fennell 
Picatinriy  Arsenal.  NJ  07806-5000 

8  Commander 

U.S.  Army  Annament  Research, 
Development,  and  Engineering  Center 
ATTN:  SMCAR-AE,  J.  Picard 
SMCAR-AEE-B, 

A.  Bcardcll 
D.  Downs 

S.  Einstein 

A.  Bracuti 
D.  Chiu 

SMCAR-AEE,  J.  Lannon 
SMCAR-AES,  S.  Kaplowitz 
Picaiinny  Arsenal,  NJ  07806-5000 

7  Commander 

U.S.  Army  Armament  Research, 
Development,  and  Engineering  Center 
A'rTN:  SMCAR-FSA-T,  M.  Salsbury 
SMCAR-FSE, 

T.  Cora 

B.  KnuteLsky 
K.  C.  Pan 
W.  Davis 

C.  Durham 
A.  Graf 

Picaiinny  Arsenal,  NJ  07806-50(X) 

6  Commtmdcr 

U.S.  Army  Armament  Research, 

Development,  and  Engineering  Center 
AOl’N:  SMCAR-EG-, 

G.  Ferdinand 

H.  Naber-Libby 
R.  Lundbc.’'g 

J.  Niles 
R.  Morcira 
W.  Moreili 

Picaiinny  Arsenal,  NJ  07806-5000 


No.  of 

Copies  Organization 

2  Commander 

U.S.  Army  Research  Office 
ATFN:  Technical  Library 

D.  Mann 
P.O.  Box  12211 

Rc.scarch  Triangle  Park,  NC  2'’709-221 1 

1  Commander 

U.S.  Army  Belvoir  R&D  Center 
ATTN:  STRBE-WC, 

Technical  Library  (Vault) 

Bldg.  315 

Fort  Belvoir,  VA  22060-5606 

1  Commander 

U.S.  Army  TRAC  -  Fort  Ixe 
Defense  Logistics  Studies 
Fort  Lee,  VA  23801-6140 

1  President 

U.S.  Army  Artillery  B,„u'd 
Fort  Sill,  OK  73503 

1  Commandant 

U.S.  Army  Command  and  General  Staff  College 
Fort  Leavenworth,  KS  66027-5200 

1  Commandant 

U.S.  Army  Special  Warfare  School 
ATl’N:  Rev  and  Tng  Lit  Div 
Fort  Bragg,  NC  28307 

1  Commander 

Radford  Army  Ammunition  Plant 
ATTN:  SMCRA-QA/Hl  Library 
Radford,  VA  24141 

1  Commandant 

U.S.  Army  Field  Artillery  School 
AinM;  STSF-TSM-CN 
Fort  Sill,  OK  /3503-5600 

4  Deputy  Commander 

Strategic  Dcfcn.se  Command 
AnX:  SFAE-SD-HVL, 

S.  Smitii 
LTC  Kce 
D.  Lianos 
T  Aden 
P.O.  Bo.\  15(X) 

Huiit.svillc,  AL  35887-8801 
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3  Commander 

U.S.  Army  Foreign  Science  and  Technology 
C  enter 

ATTN:  AMXST-MC-3. 

S.  LcBcau 
C.  Boiler 
220  Seventh  St„  NE 
Charlottesville,  VA  22901 


1  Commandant 

U.S.  Army  Field  Artillery  Center  and  School 
Arm:  ATSF-CO-MW,  B.  Willis 
Fort  Sill.  OK  73503 


1  Naval  Sea  System  Command 
Department  of  the  Navy 
ATTN:  CSEA,  CDR  Dampicr 
06KR12 

Washington,  DC  20362-5101 


3  Commander 

Naval  Surface  Warfare  Center 
Indian  Head  Division 
ATTN:  610,  C.  Smith 
6110J,  K.  Rice 
61  IOC,  S.  Peters 
Indian  Head.  MD  20640-5035 


6  Commander 

Dahlgrcn  Division 

Naval  Surface  Warfare  Center 

AITN:  CodeG33, 

T.  Doran 
J.  Copley 

Code  U30,  Guns  and  Munitions  Division 
Code  G301,  D.  Wilson 
Code  G32,  Guns  Systems  Division 
Code  E23,  Technical  Library 
Dahlgrcn,  VA  2244S-50(X) 


1  Office  of  Naval  Research 
ATTN:  Code  473,  R.  S.  Miller 
800  N.  Quincy  St. 

Arlington,  VA  22217 

2  Commander 

Naval  Sea  Systems  Command 
ATTN:  SEA  62R 
SEA  64 

Washington,  DC  20362-5101 

1  Commander 

Naval  Air  Systems  Command 
ATTN:  AIR-954 -Technical  Library 
Washington,  DC  20360 

1  Naval  Research  Labtrratory 
Technical  Library 
Washington,  DC  20375 

2  Commander 

Naval  Surface  Warfare  Center 
AH-N:  J.  P.  Con.saga 
C.  Gotzmer 

Silver  Spring,  MD  2()902-5(XX) 

2  Commander 

Naval  Surface  Warfare  Center 
ATTN:  K.  lvim/C(Hle  R-13 

R.  Bcmcckcr/Cotic  R-13 
Silver  Spring,  MD  20902-50(X) 


1  Commander 

Naval  Surface  Weapon  Center 
Indian  Head  Divisioii 
ATTN:  Code  270P1,  Mr.  Ed  Chan 
101  Straus  Avenue 
Indian  Head,  MD  20640 

1  Commander 

Naval  Surface  Weapon  Center 
Indian  Head  Division 
ATTN:  Code  3120,  Mr.  Robert  Rasi 
101  Straus  Avenue 
Indian  Hc;)d,  MD  20640 

1  Commander 

Naval  Surface  Weapon  Center 

Indian  Head  Division 

ATTN:  Code  210P1,  Mr.  Ron  Simmons 

101  Straus  Avenue 

Indian  Head,  MD  20640 

2  Commander 

Naval  Surface  Wcaprin  Center 
Indian  Head  Division 
ATTN:  Code  6210, 

Sharon  Boyles 
Norberto  Almcyda 
101  Straus  Avenue 
Indian  Head,  MD  20640 
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1  Naval  Air  Warfare  Center 

ATTN:  Code  3891,  Mr.  Chan  Price 
China  Lake,  CA  93555 

1  Commander 

Naval  Air  Warf;irc  Center 

ATTN:  Code  3891,  Ms.  Alice  Atwood 

China  Lake,  CA  93555 

1  Commander 

Naval  Weapons  Center 
Amvl:  Code  388,  C.  F.  Price 
Info  Science  Div 
China  Lake,  CA  93555-6001 

1  OL,AC  PUTSTL 

ATTN:  D.  Shiplett 
Edwards  AFB,  CA  93523-5000 

10  Central  Intelligence  Agency 

Office  of  Ccnixal  Rcl'crcnce  Di.sseiniriation 
Branch 

Room  GE-47  HQS 
Washington,  DC  20502 

1  Central  Intelligence  Agency 
ATTN:  Joseph  E.  Backofen 
HQ  Room  5F22 
Washington,  DC  20505 

5  Director 

Sandia  National  Laboratories 
ATTIS:  T.  Hitchcock 

R.  Woodfin 
D.  Bcn.son 

S.  Kempka 
R.  Beasley 

Advanced  ProjecLs  Div  14 
Organi/ation  912.5 
Albuquerque,  NM  87185 

2  Director 

Los  Alamos  National  Laboratory 
ATTN:  B.  Kaswhia 
H.  Davis 

Los  Alamos,  NM  87545 


1  Director 

Lawrence  Livermore  National  Laboratory 
ATTN:  M.  S.  L-355,  A.  Buckingham 
P.O.  Box  808 
Livermore,  CA  94550 

2  Director 

Sandia  National  Laboratories 
Combustion  Rc.search  Facility 
ATT'N:  R.  Armstrong 
S.  Vosen 
Division  8357 

Livermore,  CA  94551-0469 

1  University  of  Illinois 

Dept,  of  Much , /Indust.  Engr. 

ATTN:  Professor  Herman  Kricr,  144  MEB 
1206  N.  Green  St. 

Urbana,  IL  61801 

1  The  Johns  Hopkins  University./CPIA 
ATTN:  T.  Christian 

10630  Little  Patuxent  Parkway,  Suite  202 
Columbia,  MD  21044-3200 

2  Pennsylvania  State  University 
Dept,  of  Mechanical  Engr. 

ATTN:  Jeff  Brown 

312  Mechanical  Engineering  Bldg. 
University  Park,  PA  16802 

1  North  Carolina  State  University 
ATTN:  John  G.  Gilligan 

Box  7909 

1110  Burlington  Engineering  l^ibs 
Raleigh,  NC  27695-7909 

2  Institute  for  Advanced  Studies 
ATTN:  Dr.  H.  Fail 

Dr.  T.  Kiehne 
4030-2  West  Baker  Lane 
Au.stin,  TX  78  759-5329 

1  SRI  Intcmutinnal 

Propulsion  Sciences  Divi.sion 
ATTN:  Technical  Library 
333  Ravenswood  Avc. 

Menlo  Park,  CA  94025 
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1  SPARTA 

ATTN:  Dr.  Michael  Holland 
9455  Townc  Center  Dr. 

San  Diego,  CA  92121-1964 

5  FMC  Corporation 

ATTN:  Mr.  G.  Johnson 
Mr.  M.  Scale 
Dr,  A.  Giovanetti 
Mr.  J.  Dyvik 
Dr.  D,  Cook 
4800  East  River  Rd. 

Minneapolis,  MN  55421-1498 

2  Hercules  Inc. 

Radford  Army  Ammunition  Plant 
Manager  Manufacturing  Engineering 
Department 

ATTN:  D.  A.  Worrell 
Edward  Sanford 
P.O.  Box  1 
Radford,  V A  24141 

1  Hercules  Inc. 

ATTN:  Dr.  Richard  Cartwright 
100  Howard  Blvd. 

Kenvil,  NJ  07847 

3  GT  Devices 

ATTN:  Dr.  S.  Goldstein 
Dr,  R.  J.  Grcig 
Dr.  N.  Winsor 

5705A  General  Washington  Dr. 
Alexandria,  VA  22312 

3  General  Dynamics  Land  Systems 

ATTN:  Dr.  B.  VanDeusen 
Mr.  F.  Lunsford 
Dr.  M.  Weidner 
P.O.  Box  2074 
Warren,  MI  4S()9()-2{)74 

2  Alliant  Techsystems,  Inc. 

ATTN:  R  E.  Tompkins 

J.  Kennedy 
7225  Northland  Dtive 
Brooklyn  Park,  MN  55428 


4  01  in  Ordnance 

ATTN:  V.  McDonald,  Librai7 
Hugh  McElroy 
Mr.  Thomas  Bourgeois 
Mr.  Dennis  Wortfiington 
P.O.  Box  222 
St.  Marks,  FL  32355 

1  Paul  Gough  A.s.sociatcs,  Inc. 

ATTN;  P.  S.  Gough 
1048  South  St. 

Portsmouth,  NH  03801-5423 

1  Physics  International  Library 
ATTN:  H.  Wayne  Wampler 
P.O.  Box  5010 

San  Leandro,  CA  94577-0599 

2  Rockwell  intemalional 
Rocketdyne  Division 
ATTN:  BA08, 

J.  E.  Flanagan 
J.  Gray 
6633  Canoga  Avc. 

Canoga  Park,  CA  91304 

2  Princeton  Combustion  Research 

Laboratories,  Inc. 

ATTN:  M.  Summcrficld 
N.  Messina 

Princeton  Corporate  Plaza 
11  Dccipark  Drive 
Bldg.  IV,  Suite  1 19 
Monmouth  Junction,  NJ  08852 

2  Science  Applications,  Inc. 

ATTN;  J.  Batteh 

L.  Thornhill 
1519  John.son  Ferry  Rd. 

Suite  3(K) 

Marietta,  GA  3(XX)2-(>438 

1  Eli  Freedman  <&  Associates 

ATTN:  E.  Freedman 
2411  Diana  Rd. 

Baltimore,  MD  21209 
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1  Rockctdync 

ATTN:  Mr.  Oito  Hcincy 
Mail  Stop  BA26 
6633  Canoga  Avenue 
Conoga  Park,  CA  913(M 

1  Thiokol  (Longhorn  Division) 

ATTN:  Dr.  David  Dillchay 
Mail  Stop  703-11 
P.O.  Box  1 149 
Marshall,  TX  75671 

1  Thiokol  (Elkton  Division) 

ATTN:  Dr.  Rodney  Wilier 
55  Thiokol  Road 
Elkton,  MD  21922 

1  Veritay  Technology,  Inc. 

ATTN:  Mr.  E.  Fisher 
4845  Millcr.spoit  Highway 
Ea-st  Amherst,  NY  14051-0305 

1  Veritay  Technology,  Inc. 

4845  Millcrsport  Hwy. 

P.O.  Box  305 

East  Amherst,  NY  14051  •0305 

1  Battel  Ic 

ATTN:  TACTEC  Library,  J.  N  Huggins 
505  King  Avc. 

Columbus,  UH  43201-2693 

2  California  Institute  of  Technology 
Jet  Propulsion  Lalxiratory 
ATTN:  L.  D.  Strand,  MS  125-224 

D,  Elliot 

48(X)  Oak  Grove  Dr. 

Pa.Siidcna,  CA  91 109 

1  General  Electric  Co. 

Dcfcn.se  Systems  Division 
AITN:  Dr.  J.  Mand/.y 
Mail  Drop  43-220 

100  Pla.stics  Avc. 

Pittsfield,  MA  01201 

2  SAIC 

ATTN:  Mr.  N.  Sinha 
Dr.  S.  Dash 

501  Office  Center  Drive 

Fort  Washington,  PA  19034-3211 


No.  of 
Comes 

1 


1 


1 


1 


1 


4 


Organi/ation 

State  University  of  New  York 
Department  of  Elecmical  Engineering 
ATTN:  Dr.  W.  J.  Sargeant 
Bonner  Hall  -  Room  312 
Buffalo,  NY  14260 

Science  Application  Intcmatiimal  Coqioration 
ATTN:  Dr.  George  Chry.s.simicll is 
8400  Normandclc  Blvd. 

Suite  939 

Minncapttlis,  MN  55437 
Hercules,  Inc. 

ATTN:  Dr.  Richard  Cartwright 
100  Howard  Blvd. 

Kenvil,  New  Jersey  078'',7 

Sandia  National  Laboratories 
ATTN;  Mr,  Mark.  Grulxilich,  DIV  2515 
P.O.  Box  58(X) 

Albuquerque,  NM  87185 

IMl  Services  USA 
ATTN:  Mr.  G.  Rashba 
2  Wisconsin  Circle 
Suite  420 

Chevy  Cha.se,  MD  20815 


Aberdeen  Proving  Ground 

Cdr,  USACSTA 
ATIN:  S.  Walton 
G.  Rice 
D.  Lacey 
C.  Hcrud 
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2  RARDE 

GS2  Division 
Building  R31 
ATTN:  Dr.  C.  Woodley 
Dr.  G.  Cook 
Fort  Halstead 

Sevenoaks,  Kent  TNI 4  7BP 
ENGLAND 

1  Materials  Research  Laboratory 

Salisbury  Brunch 

AITN:  Anna  Wildcgger  Gaissmaicr 
Explosives  Ordnance  Division 
Salisbury  South  Australia  5108 

1  Laboraturio  Quiniico  Ccnual  dc  Arniamcnto 

ATTN;  Captain  Juan  F.  Hernandez  Tamayo 
Apartado  1105 
28080  Madrid 
SPAIN 

1  R&D  DepajttneiU 

ATTN:  Dr.  PieiTC  Archambaull 

5  Montec  dcs  Arsenaux 

Le  Gardcur,  Quebec,  Canada  J5Z  2P4 

1  TZN  Forschungs-  uiid  Entwicklungszcnuuin 

A1TN:  Mr.  T.  Weise,  Dr.-Ing. 

Abt.  SE,  Fachbcreich; 
Hochleistungspulstixhnik 
Neuensothriether  Stralle  20 
D-3104  Untcrlilb,  Germany 

1  Emst-Mach-Institut 

A'lTN:  Dr.  Gustav-Adolf  SchrOdcr 
HaupisUalk;  18 

D-7858  Weil  am  Rhein,  Germany 

2  Institui  Hrai'Co  Allcmand 
ATTN:  Dr.  M.  Sumiranl 

Mr.  D.  Grune 

F  68301  SAINT-LOUIS  Cddcx,  12, 
rue  de  ITndu.stric,  B.P.  301,  France 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
j  comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

j  1.  ARL  Report  Number  ARL-MR-6 1 _ Date  of  Report  April  1993 _ 

2.  Date  Report  Received _ 

■■  i 

I  3.  Docs  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Infomation  .source,  design  data,  procedure,  source  of 
ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  a.s  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 

'  I  _ _ _ 

r  j  6.  General  Comments.  .  What  do  you  think  should  be  changed  to  improve  futuic  reports?  (Indicate 

.  ]  changes  to  organization,  technical  content,  fonnat,  etc.) _ 


Organization 


CURRENT  Name 

ADDRE.SS  _ 

Sheet  or  P.0,  Box  No. 


City,  State,  Zip  Code  j 

7.  If  indicating  a  Change  of  Addrcs.s  or  Address  Correction,  please  provide  the  Current  or  Correct  addic.ss 
above  and  the  Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

(Remove  lliis  slicct,  fold  as  indicated,  tajx;  closed,  and  mail.) 
(DO  NOT  STAPLE) 


